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ABSTRACT 

Aims. XO-2 is the first confirmed wide stellar binary system where the almost twin components XO-2N and XO-2S have planets, 
and it is a peculiar laboratory to investigate the diversity of planetary systems. This stimulated a detailed characterization study of the 
stellar and planetary components based on new observations. 

Methods. We collected high-resolution spectra with the HARPS-N spectrograph and multi-band light curves. Spectral analysis led to 
an accurate determination of the stellar atmospheric parameters and characterization of the stellar activity, and high-precision radial 
velocities of XO-2N were measured. We collected fourteen transit light curves of XO-2Nb used to improve the transit parameters. 
Photometry provided accurate magnitude differences between the stars and a measure of their rotation periods. 

Results. The iron abundance of XO-2N was found +0.054 dex greater, within more than 3<x, than that of XO-2S. The existence of a 
long-term variation in the radial velocities of XO-2N is confirmed, and we detected a turn-over with respect to previous measurements. 
We suggest the presence of a second massive companion in an outer orbit or the stellar activity cycle as possible causes of the observed 
acceleration. The latter explanation seems more plausible with the present dataset. We obtained an accurate value of the projected 
spin-orbit angle for the XO-2N system (/1=7°±11°), and estimated the real 3-D spin-orbit angle (i/r=27degrees). We measured the 
XO-2 rotation periods, and found a value of P=41,6± 1.1 days in the case of XO-2N, in excellent agreement with the predictions. The 
period of XO-2S appears shorter, with an ambiguity between 26 and 34.5 days that we cannot solve with the present dataset alone. 
The analysis of the stellar activity shows that XO-2N appears to be more active than the companion, and this could be due to the 
fact that we sampled different phases of their activity cycle, or to an interaction between XO-2N and its hot Jupiter that we could not 
confirm. 

Key words. (Stars:) individual: XO-2S, XO-2N — Stars: fundamental parameters, abundances — planetary systems — techniques: 
radial velocities, photometric 
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1. Introduction 

The diversity of the more than 1,700 extrasolar planets dis¬ 
covered so fafl (~850 of which are in multi-planet systems; 
iRowe et al.ll2014l) . from their orbital architectures to the astro- 
physical environments where they reside, represents a very com¬ 
plex issue to be investigated to understand the mechanisms of 
their formation and dynamical evolution. 

In the era of comparative exo-planetology, a fundamental ap¬ 
proach for the understanding of the different exoplanetary prop¬ 
erties relies on the characterization of planet-host stars and statis¬ 
tical analysis of the planetary system frequency. Several studies 
shown that the properties of the planetary system architectures 
and physical characteristics of the planets depend upon stellar 


prope rties, such as the mass dJohnson_etaL| 


2011 


201 3j) and the metallicit y dSozzetti et alJ 


201C; iBonfih^eyd, 


2009; 


Santos et al 


iMortier et al.ll2012h . Moreover, during the last years the 
search for exoplanet s was extended in envir onments like giant 
stars (e.g. Kepler-91. lLillo-Box et al.ll2014bllah .open and globu¬ 
lar cluster s, and only recent l y planetary compani ons were found 
(see, e.g.,[Ouinn et al. 2012; Meib om et alJl2013l and references 
therein). 

Binary systems are interesting astrophysical environments to 
search for planets. Considering that nearly half of the solar-type 
stars are gra vitationally bound with at least another stellar com¬ 
panion (e.g. lRaghavan et aPl2010l : lDuchene & Kraus|l2013l) . bi¬ 
nary systems naturally represent a typical environment to be ex¬ 
plored for studying the processes of planet formation and evolu¬ 
tion leading to very different planetary architectures. 

Planet formation in a binaries has been considered by var¬ 
ious authors. For coplanar orbits be tween the disk and the 


binary the general expectation (e.g ., Marzari & Scholll [2000; 


iKlev & Nelsonl l2008t iMarzari et al J l2(B 2l) is that the tidal ef¬ 
fect of the companion would induce an eccentricity growth in 
the planetesimal population that might inhibit planet formation 
if the binary separation is < 50 AU. In this configuration even 
the formation of icy planetesimals is jeopardized by the high 
temperatures in the discs and the form ation of breaking waves 
dNelsonll2000t iPicogna & Marzarill201 3l) . The situation is more 
complicated for misaligned systems where the planetesimals and 
gas discs precess around the binary orbital plane. Planetesi¬ 
mals naturally develop strong differential precession that would 
inhibit planet formation by increasing the velocity dispersion 
dMarzari et al.112009 1. On the other hand, if the gas disk is ra¬ 
dially narrow (which is expected for separation of the order of 
50 AU because of tidal truncation), it will precess rigidly and 
would drag the planetesimalsjnto rigid precession, mitigating 
the previous effect dFragner et al.ll201lh . unless the inclination is 
so large for Lidov-Kozai effects to become important. Also in 
this case, the effects should be relevant mostly for binary separa¬ 
tions smaller than 50-100 AU. For wider, misaligned binaries the 
gas disk on the other hand is expected to reach a_ quasi-steady, 
non precessing warped configuration tsee iFacchini et all [20131 


* Based on observations made (/) with the Italian Telescopio 
Nazionale Galileo (TNG), operated on the island of La Palma by the 
INAF - Fundacion Galileo Galilei (Spanish Observatory of Roque de 
los Muchachos of the IAC); («) with the Copernico and Schmidt tele¬ 
scopes (INAF - Osservatorio Astrofisico di Padova, Asiago, Italy); (Hi) 
with the IAC-80 telescope at the Teide Observatory (Instituto de As- 
trofisica de Canarias, IAC); (tv) at the Serra la Nave "M.G. Fracastoro" 
Astronomical Observatory (INAF - Osservatorio Astrofisico di Cata¬ 
nia); (v) at the Astronomical Observatory of the Autonomous Region of 
the Aosta Valley (OAVdA). 

1 NASA exoplanet archive, http://exoplanetarchive.ipac.caltech.edu/ 


for the case of a circumbinary disk). If the warp is located in 
the planet formation region, it might affect the process, although 
such effects have not yet been studied. 

Several examples of planets orbiting only one of 
the binary components, defin e d as S-type planets , are 

2012 ; iMuerauer et al.l 12014 and 
2014 for a recent review). The 


known (e.g. [Roell et al 


IThebault & Haghighipour 


existence of S-type planets rises several questions: how the 
presence of a stellar companion can affect the formation, 
survival and dynamical evolution of such planets? Is it possible 
to identify which properties of a stellar system most influence 
the physical characteristics of the hosted planets? Does any 
constraint exist with respect to the case of isolated stars, or 
any dependence from the physical separation of the binary 
components? Observational difficulties e xist in searching for 
S-typ e planets in multiple stellar systems (lEggenberger & Udrvl 
I2010I) . in particular those with a very small angular separation 
(less then ~2"), since, for instance, the components cannot be 
observed as two isolated targets in the spectrograph’s fiber or 
slit. In many cases of confirmed extrasolar planets the presence 
of a stellar companion to the host_star was discovered after the 
detection of the planet (e.g. [Roell et alJl2012 ). Binarity poses a 
major challenge in understanding how such planets could exist 
and their evolutionary history in the stellar system. 

Very wide binaries, with semi-major axes of the order of 10 2 3 
AU and sky projected angular separations of several arcsec, of¬ 
fer the best opportunity and less observational complications to 
search for planets orbiting each component. At the same time 
they allow an accurate determination of the properties of the in¬ 
dividual stars. Despite the large separation between the stellar 
components, these systems are an interesting subject for dynam¬ 
ical studies also related to the presence of planets. 

iKaib et all (1201 3l) and iKaib & Ravmondl (120141) studied the 
perturbations produced by passing stars and the tidal effects of 
the Milky Way on the orbits of very wide binary stars. Interest¬ 
ingly they found that occasionally the orbits could become ex¬ 
tremely eccentric, resulting in the collision of the binary com¬ 
ponents or, less dramatically, they could evolve toward close 
or contact systems. Moreover, simulations bv lKaib et al. ( 2013 ) 
indicate that very wide binary companions may often strongly 
reshape planetary systems after their formation, by causing the 
ejection of planets from the system or heavily changing the or¬ 
bital eccentricities of those that survive. Another intriguing ex¬ 
ample of a plausible fate for planets in wide binarie s is th e 
’bouncing’ scenario investigated by iMoeckel & Veras! ( 2012 ). 
They found that in 45 to 75 per cent of the cases a planet, ini¬ 
tially orbiting one of the binary components and then ejected 
from its native system by planet-planet scattering, could jump 
up and down passing from the gravitational influence of one star 
to the other, for more than a Myr. This situation may trigger or¬ 
bital instability among existing planets around the companion, 
and in some cases result in an exchange of planets between the 
two stars, when both host multiple planetary systems. 

At present, few very wide binaries are known where a plan¬ 
etary system was discovered around one of the stellar compo¬ 
nents. The star HD20782 is the binary companion of HD20781, 
with a projected separation of -9000 AU, and it hosts a Jupiter- 
mass planet on a very eccentric orbit at ~1.4 AU (1 Jone s et alJ 
l2006lrl . The very high eccentricity reported in the literature (e 
= 0.97) could be the result of significant perturbations experi- 


2 HD20781 should_host two Neptune-mass planets within ~0.3 AU, 

as reported by iMavor et all d201 if) , but the discovery has not yet been 
confirmed and published. 
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enced by the planet during the evolution of the system and pos¬ 
sibly due to the stellar companion. In this binary system, where 
both the components can be well analysed separately, a compar¬ 
ative study of the physical properties of the pair could be of help 
for explaining the nature of the planetary system. To investigate 
any existing link between the planet formation mechanisms and 
the chemical composition of the host stars, iMack et aD (2014 ) 
performed a detailed elemental abundance analysis of the atmo¬ 
spheres of HD20781/82, by considering both stars to host plan¬ 
ets. This kind of study, when at least one component hosts a 
close-in giant planet, could result in the evidence of chemical 
imprints left in the parent star suggesting the ingestion of mate¬ 
rial from the circumstellar disk (and possibly also of planetary 
origin) driven by the dynamical evolution of the planet orbit. In 
fact, when a star with a close-in giant planet is found to be en¬ 
riched with elements of high condensation temperature, as sug¬ 
gested by lSchuler et al.l (1201 lbi) . this can be related to the inward 
migration of the planet from the outer regions of the circumstel¬ 
lar di sk, where it formed, to the present p o sition closer to the 
star (llda & Linll2008b iRavmond et alJl20llh . IMack et all (1201 4l) 
found for both stars a quite significant positive trends with the 
condensation temperature among the elemental abundances, and 
suggest that the host stars accreted rocky bodies with mass be¬ 
tween 10 and 20 M© initially formed interior to the location of 
the detected planets. Searches for abundance anomalies possibly 
caused by the ingestion of planetary material by the central star 
were also conducted bv lDesidera et al.l (20041120071) (and refer¬ 
ences therein) for a sample of wide binaries. These authors found 
that the amount of iron accreted by the nominally metal richer 
companion (in binaries with components having T e ff>5500 K) 
is comparable to the estimates of the rocky material accreted by 
the Sun during its main-sequence lifetime, and therefore con¬ 
cluded that the metal enrichment due to the ingestion of material 
of planetary origin should not be a common event. 

A representative case of a very wide binary with one com¬ 
ponent hosting a close-in giant planet is that of the XO-2 system 
(projected separation ~31"), where the star XO-2N L orbited by 
the transiting hot Jupiter XO-2Nb ( iBurke etakl 1 2007 1. An ele¬ 
mental abundance analyses of XO-2N and its companion XO-2S 
was performed bv lTeske et al.l ( 2013 ). who determined the stel¬ 
lar abundances of carbon and oxygen, as well as iron and nickel. 
Their goal was to probe the potential effects that planet forma¬ 
tion and evolution might have had on the chemical composition 
of XO-2N, where the companion XO-2S was treated as a non¬ 
hosting planet star. Later, iDesidera et al.l (120141) discovered that 
the star XO-2S actually hosts a planetary system formed by a 
planet slightly more massive than Jupiter orbiting at 0.48 AU 
and a Saturn-mass planet at 0.13 AU (with cautious evidence of 
a long-term trend in the radial velocities time series possibly due 
to an outer companion, yet of unknown nature). The discovery is 
of particular relevance because it represents the first confirmed 
case of a very wide binary in which both components host plan¬ 
ets (see Fig. [I] for a sketch of the XO-2 planetary systems). The 
XO-2 binary offers a unique opportunity to explore the diver¬ 
sity of planet formation mechanisms, by considering that, while 
the parent stars are almost equal in their main physical proper¬ 
ties, the planetary systems are different. This finding has conse¬ 
quently motivated a detailed comparative study of the whole sys¬ 
tem, which is the subject of this paper, through the use of new 
high-resolution and high S/N HARPS-N spectra and dedicated 
differential photometry. It is crucial to highlight any existing dif¬ 
ference between the properties of the two stars which could have 
played a distinctive role in the formation processes of the plan¬ 
etary systems: any observed difference, if correctly interpreted. 


can tell a part of that story. While the analysis of the XO-2 sys¬ 
tem was still in progress, the discovery of a hot Jupiter around 
each of the twins, metal-rich ([Fe/H ]~0.25 dex), F-type compo¬ 
nents of the WASP-94 wide binary (projected separation ~2700 
AU) was announced dNeveu-VanMalle et al.ll2014b . This system, 
with host star properties and planetary system architectures dif¬ 
ferent from those of XO-2, indeed makes more intriguing the 
understanding of the role of stellar multiplicity on the planet for¬ 
mation and evolution. 

This paper is organized as follows. We first describe in Sect|2] 
the spectroscopic and photometric datasets used in this study. We 
present in Sectj3]a new determination of the basic stellar param¬ 
eters based on HARPS-N spectra, as well as results of a detailed 
and homogeneous analysis aimed at searching for differences 
in the iron abundance between the XO-2 components with high 
level of confidence. In Sectj4]we present and discuss results from 
a photometric follow-up of the two stars, including the measure¬ 
ment of the stellar rotation periods and the modelling of a new 
dataset of 14 transits of the planet XO-2Nb. Sectj3] is dedicated 
to the XO-2N planetary system. There we present a new analysis 
of the radial velocity time series, thanks to which we confirm the 
existence of an acceleration with a not yet established cause. In 
the same Section we also present new observation and analysis 
of the Rossiter-McLaughlin effect, better constraining the value 
of the projected spin-orbit angle A. A detailed analysis of the stel¬ 
lar activity for the XO-2 components is presented and discussed 
in Sectj6] We finally conclude with Sect[7] where we speculate 
about the evolution and long-term stability of the XO-2 planetary 
systems. 

2. Observations and data reduction methods 

2.1. Spectroscopy 

The spectra of the XO-2 components analysed in this work 
were collected with the high-resolution HARPS-N spectrograph 
(iCosentino et al.l 12012 ) installed at the Telescopio Nazionale 
Galileo (TNG) on La Palma (Canary islands). The observa¬ 
tions were carried out in the framework of the larg e programme 
Global Architecture of Pla netary Systems fGAPS : ICo vino et aD 
l20H IDesidera et al.ll20H . 

For the star XO-2N we collected 43 spectra between Novem¬ 
ber 20, 2012 and October 4, 2014, while the companion XO-2S 
was observed at 63 individual epochs between April 21, 2013 
and May 10, 2014. The Th-Ar simultaneous calibration was not 
used to avoid contamination by the lamp lines. This has no sig¬ 
nificant impact on the measurements of the radial velocity, be¬ 
cause the drift correction with respect to the reference calibra¬ 
tion shows a dispersion of 0.8 m s -1 , which is smaller than the 
median radial velocity (RV) uncertainty due to photon-noise, i.e. 
2.0 m s _1 for XO-2N and 2.2 m s _1 for XO-2S. The reduction 
of the spectra and the RV measurements were obtained for both 
components using the latest version (Nov. 2013) of the HARPS- 
N instrument Data Reduction Software pipeline and applying a 
K5 mask. The measurement of the RVs is based on the weighted 
cross-correlatio n function (CCF) method (Bara nne et al.lll996t 
iPene et akll2002l) . 

2.2. Photometry 

We collected and analysed new photometric light curves for both 
components of the XO-2 system with three different facilities. In 
particular, XO-2S was intensively monitored with the HARPS- 
N spectrograph. When it became clear that the star had plane- 
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Fig. 1 . Comparison between the planetary systems orbiting the XO-2 
binary stars and the architecture of the inner Solar System. The case of 
XO-2N is shown in the upper panel, while that of XO-2S is represented 
in the lower one. In both plots, the red ellipses (solid lines) represent the 
orbits of the XO-2 planets. The orbits of Mercury, Venus and the Earth 
are indicated with black, blue and green dotted lines respectively, also 
labelled with the planet names. The cross symbol in the middle of each 
plot indicates the star location. 


2.2.1. The APACHE dataset 

A follow-up of the XO-2 field was conducted for 42 nights be¬ 
tween December 2, 2013 and April 8, 2014 with one of the 40cm 
telescopes composing the APACHE array dSozzetti et alJl2013h . 
based at the Astronomical Observatory of the Autonomous Re¬ 
gion of the Aosta Valley (OAVdA, +45.7895 N, +7.478 E, 1650 
m a.s.l.). Each telescope is a Carbon Truss f/8.4 Ritchey-Chrtien 
equipped with a GM2000 10-MICRON mount and a FLI Pro¬ 
line PL1001E-2 CCD Camera, with a pixel scale of 1.5 "/pixel 
and a field of view of 26'x26'. The observations were carried 
out using a Johnson-Cousins I filter. The images were reduced 
with the standard pipeline TEEPEE written in IDL0 and regu¬ 
larly used for the APACHE purposes (see Giacobbe et al.ll2012h . 
The APACHE data were used to tentatively derive the rotation 
period of the two stars and to monitor their magnitude difference 
in I band, as described in Sect|4] 

2.2.2. The TASTE dataset 

Since February 2011 several transits of the hot Jupiter XO-2Nb 
have been observed at high temporal cadence with the 182cm 
and the Schmidt 92/67 telescopes at the Asiago Astrophysical 
Observatory (+45.8433 N, +11.5708 E, 1366m a.s.l.), in the 
framework of the TASTE Proj ect (The Asiago Search for Transit 
time variations of Exonlanets: lNascimbeni et al.M201lb . A subset 
of these measurements was used to estimate the magnitude dif¬ 
ferences of the two XO-2 binary components in the R, band. 
Two additional transits were also collected with the IAC-80 tele¬ 
scope at t he Teide Observatory, u sing the CAMELOT CCD im¬ 
ager (see lNascimbeni et al.ll2013l for a description of the set-up), 
and one supplementary light curve was obtained with the same 
telescope but using the T CP camera (Tromso CCD Photometer; 
lOstensen & SolheirrfoOOOt) . All these observations, for a total of 
14 light curves, were carried out in the Cousin/Bessell R band 
and were analysed together to improve the determination of the 
parameters of XO-2Nb. 

2.2.3. The Serra la Nave dataset 

The XO-2 system was also observed for five nights, between 
May 6, 2013 and May 7, 2014 with the robotic 80cm f/8 Ritchey- 
Chretien telescope APT2, operated by the INAF-Catania As¬ 
trophysical Observatory and located at Serra la Nave (SLN, 
+37.692 N, +14.973 E, 1725m a.s.l.). The telescope is equipped 
with a CCD detector Apogee Alta U9000 (3Kx3K pixels) op¬ 
erated in binning 2x2, corresponding to a pixel scale 0.76 arc- 
sec/pixel. The data were reduced with the I RAF packageQ us¬ 
ing the standard procedure for overscan, bias and dark subtrac¬ 
tion and flat fielding. The aperture photometry w as done with 
the software SExtractor (iBertin & Arnoutsll 19961) . The data of 
SLN were collected in the BVRI bands and used to measure the 
magnitude differences between the two stars. 


3. Stellar parameters 

3.1. Spectroscopic analysis of the individual stars 


tary companions, we started a dedicated photometric follow-up 
of the target to estimate the stellar rotational period, further char¬ 
acterizing its level of activity and look for possible transits of its 
exoplanets. 


The components of the XO-2 binary system have a separation 
of ~31 arcsec (corresponding to a projected distance of ~4600 
AU assuming a distance of ~ 150 pc, as estimated bv lBurke et al.l 

3 Registered trademark of Exelis Visual Information Solutions. 

4 http://iraf.noao.edu/ 
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l2007 i). and share, within the uncertainties, the proper-motion 
vector as listed in the UCAC4 catalogue. Essential information 
about the XO-2 system is presented in Table[T| A standard analy¬ 
sis of the HARPS-N spectra was performed with different meth¬ 
ods to derive the stellar atmospheric parameters of the two com¬ 
ponents, with those for XO-2S first presented in iDesidera et alJ 
(20141) . Table [2] summarizes our results, which represent the 
weighted averages of the individual measurements. We used im¬ 
plementations of both the equi valent width and the spectral syn¬ 
thesis meth ods, as describe d in Biazzo et alJ (l2012l) . ISantos et all 
(120131) . and lEsposito et al.l (I2014I) . It is interesting to note that 
the effective temperature difference between the two stars de¬ 
rived spectroscopically is in good agreement with the result ob¬ 
tained from photometric measurements. In fact, following the 
same analysis performed bv lMunari et al.l d2014l) for half-million 
of stars of the RAVE spectroscopic survey, if we assume a color 
excess E(B - V)=0.02 mag (as explained in Sect l3.3l > the XO- 
2S star results 60+33 K hotter than the companion. The stellar 
mass, radius and age were determined by comparing our mea¬ 
sured effective temperature, iron abundance, and surface gravity 
with t he Yonsei-Yale (Y-Y) evolution ary trac ks (Demarque et al.1 
2004 ) through the ^-square statistics dSanterne et al. 2011 ). The 
adopted errors include an extra 5% in mass and 3% in radius 
added in quadrature to the formal errors to take systematic uncer¬ 
tainties in stellar models into account dSouthworthl201 ll) . Taking 
advantage of the fact that the planet XO-2Nb transits its parent 
star, we derived a second, slightly more accurate estimate for the 
mass, radius and age of XO-2N, using the stellar density infor¬ 
mation as determined from the analysis of the transit light curves 
(e.g.. fSozzetti et ahll2007i) . 


3.2. Differential spectral analysis 

Thanks to the high-quality of the HARPS-N spectra, we per¬ 
formed a detailed and homogeneous spectral analysis of XO-2N 
and XO-2S to search for possible differences in the iron abun¬ 
dance between the components. The average values for [Fe/H] 
obtained from the analysis of the individual stars (Table [2} in¬ 
dicate that the North component is characterised by higher iron 
abundance, but the two estimates are compatible within the un¬ 
certainties. In order to validate this difference with a high level 
of confidence, we determined the difference in [Fe/H] very ac¬ 
curately with a dedicated analysis. Together with the spectra of 
both targets, we also acquired three solar spectra through ob¬ 
servations of the asteroid Vesta. This allowed us to perform a 
differential analysis for each component with respect to the Sun, 
thus avoiding the contribution due to the uncertainties in atomic 
parameters, such as the transition probabilities. All HARPS-N 
spectra were shifted in wavelength and then co-added in order to 
obtain very high S/N spectra (~ 180 - 250 at A ~ 6700 A). We 
performed the spectral analysis following three main steps: 

i. Analysis based on equivalent widths. We measured the 
equivalent widths (EWs) of iron lines on one-dimensional spec¬ 
tra using the splot task in IRAF, paying attention to trace as 
much as possible the same position in the continuum level for 
the spectra of both components and the asteroid Vesta. We con¬ 
sidered the line-list from iBiazzo et al . (1201 2). and adopted the 
same procedure as these authors for the differential analysis rel¬ 
ative to the Sun, the rejection criteria of bad lines, the line broad¬ 
ening mechanisms, the model atmospheres, the measurement of 
stellar parameters and iron abundances, and their corresponding 
uncertainties. We refer to that paper and to [Biazzo et alJ (1201 11) 
for detailed descriptions of the method. In brief, we used the 


version 2013 of the MOOG code ( Sneden 1973 ) and consid¬ 
ered the abfind driver, assuming local thermodynamic equilib¬ 
rium (LTE). Initial stellar parameters were set to the solar values 
(T e ff0=577O K, log go = 4.44, and £ 0 = 1-10 km/s) for both 
components. Then, the final effective temperature (T e ff) was de¬ 
termined by imposing the condition that the abundance from the 
Fe 1 lines (log n {Fe 1 )) was independent on the line excitation po¬ 
tentials; the final microturbulence it;) by minimizing the slope of 
logn(Fei) versus the reduced equivalent width ( EW/A ); and the 
final surface gravity (log g) by imposing ionization equilibrium 
(i.e. logn(Fei)=log»(Fen)). Plots of iron abundance (logn(Fe)) 
versus excitation potential (x) and reduced EW for both com¬ 
ponents are shown in Fig. [2] which shows that the correlations 
are close to zero, as required. Using the same procedure, we ob¬ 
tained logn(Fei) 0 = logn(Fen) 0 = 7.53 + 0.05 for the Sun. We 
point out that the exact values of the solar parameters are not 
crucial as we are performing a differential study for both com¬ 
ponents with respect to the Sun. Table [2 lists the final results of 
this procedure, where the errors in iron abundance include uncer¬ 
tainties in stellar parameters and in EW measurements. Results 
from this first homogeneous analysis seems to be in accordance 
with those obtained with different methods of spectral analysis 
(Tabled. 

ii. Line-by-line analysis. Since we are interested in a homo¬ 
geneous analysis, we applied the same iterative method as in 
the first step, but here we considered a strict line-by-line anal¬ 
ysis, i.e. the same lines were used for both components, with¬ 
out any implementation of line rejection criteria. This allowed 
us to avoid uncertainties related to both atomic parameters and 
measurements of equivalent widths, due to, e.g., blends and con¬ 
tinuum level. In the end, the analysis was based on 75 Fei and 
10 Fen lines. Thanks to such analysis, the errors in T e ff were 
reduced substantially, as reported in Tabled 

Hi. Differential analysis. Here, as stellar parameters of XO- 
2N, we considered the values derived in the line-by-line anal¬ 
ysis, while XO-2S was analysed differentially with respect to 
the companion, strictly using the same line set and EWs as in 
the previous step. We applied the di f ferenti al abu ndance method 
widely des cribed in iGratton et al] (12001 1) and IDesi dera et al. 
( 2004112006I) . which yields very accurate results when, in partic¬ 
ular, the two components are very similar in stellar parameters 
(as for the XO-2 components). In summary, the most important 
ingredient in this analysis is the difference in temperature be¬ 
tween the components. In our case, we take advantage of the 
knowledge of the luminosity difference between the targets, as 
we can assume they are at the same distance from us. We con¬ 
sidered as magnitude differences in different bands the mean val¬ 
ues reported in Sect l4.1l and as bolometric corrections those de¬ 
rived by using the code provided bv lCasagrande & VandenBeril 
(2014), which allows to estimate the bolometric corrections in 
several photometric bands using T e ff, [Fe/H], and logg as pri¬ 
mary inputs. Then, since the stars are on the main sequence, we 
could obtain accurate estimates of the difference A log g between 
their surface gravities (see Eq. 6 in lDesidera et al.l20061 where as 
stellar mass difference we assumed 0.01 M 0 , as given in Table 
0. This allowed us to accurately derive their temperature dif¬ 
ference (A7' e [i ) using the equilibrium of the iron ionization, be¬ 
cause the difference between the abundances provided by the Fe 1 
and Fen lines is strongly sensitive to T e ff (about 0.001 dex/K; 
see IGratton et al.i2001l) . The final atmospheric parameters of the 
secondary were derived as differences with respect to those of 
the primary by means of an iterative procedure (because log g, 
and [Fe/H] depend on the assumed temperatures). We inter¬ 
rupted the procedure considering reached the convergence when 
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the grid step-size was of 1 K in T e fj, 0.001 dex in log g, and 
0.01 km/s in No other cr clipping was implemented here. The 
final adopted differences in stellar parameters are listed in Ta- 
ble[3] while Fig.[3]shows A[Fe/H] as a function of the excitation 
potential. The under-abundance in iron and the higher effective 
temperature of XO-2S when compared to XO-2N are confirmed 
by this procedure. All differential parameters we obtained are 
affected by both internal and systematic errors. Internal error on 
AT e ff was derived from the line-by-line scatter of Fe i and Fe n 
and the errors in the other atmospheric parameters. Internal error 
in A log g includes the uncertainties in the differences in mass, 
magnitude, bolometric correction, and effective temperature. Er¬ 
ror in the microturbulence difference was obtained by summing 
in quadrature the error contributions in /■' of each component. The 


a lower iron content. To explain this deficiency, iRamfrez et al.l 
(1201 ll) suggest that an early depletion of metals happened dur¬ 
ing the formation phase of the 16 Cyg Bb planet. An alterna¬ 
tive and suggestive possibility is that 16 Cyg A may have in¬ 
gested a massive planet that enriched the star with iron, while 
the large orbit of 16 Cyg Bb (with semi-major axis of 1.68 AU) 
likely prevented any mass loss from the planet. We note that, 
despite A[Fe/H] is similar with that we found for XO-2, the 
amount of iron involved should be different. In fact, the XO-2 
stars are cooler than 16 Cygni A and B and have therefore more 
massive convective envelopes, implying that more iron should 
be necessary to pollute the XO -2N photosphere and prod uce al¬ 
most the same A[Fe/H]. Also lLaws & Gonzalez d200 ll) found 
the primary of the 16-Cyg system enhanced in Fe relative to 


account the uncertainties in all stellar parameters added quadrat- 

ing planets (see, e.g., Gratton et al. 2001; Desidera et al 

2004 

ically. All errors in the differential analysis are very low, as sum- 

2006; Schuler et al. 201 Teske et al. 20131 Liu et al. 

2014 


marized in Table [3 and this demonstrates how this procedure 
is efficient to unveil small differences in atmospheric parame¬ 
ters, removing many sources of systematic errors. In fact, as re¬ 
marked bv lDesidera et al.l d2006l) . external uncertainties, due to, 
e.g., wrong estimations of parallax and mass, or inadequacies of 
model atmospheres and LTE deviations, have negligible effects 
on our results, because of the similar characteristics of the com¬ 
ponents. In the end, the difference in iron abundance of XO-2N 
with respect to XO-2S is of +0.054 dex at a more than 3<x level. 

The difference in iron abundance between the two XO-2 stel¬ 
lar components poses an interesting issue. They belong to a vi¬ 
sual binary and, as normally assumed for such systems, they 
should share the same origin and initial bulk metallicity. A rel¬ 
evant characteristic of this system is that both of the stars host 
planets. Thus, for components of wide binaries where at least 
one star has a planet, a reasonable hypothesis to explain any 
measured and significant difference in their present-day elemen¬ 
tal abundances is that most likely the planet formation process 
had played a relevant role. The higher iron abundance of XO- 
2N when compared to XO-2S might be due to past ingestion 
of dust-rich or rocky material, coming from the inner part of 
the proto-planetary disk and pushed into the host star by the hot 
lupiter XO-2N as it migrated inward to its current orbit. A sec¬ 
ond mechanism, which acts on a different time scale and after the 
pre-main-sequence stage, is presented bv lFabrvckv & Tremaine! 
(20071) . who discuss the pollution of the stellar photosphere with 
metals produced by mass loss of inwar d migr ating hot Jupiters. 
Also, with their simulations iKaib et al.l ( 2013 1 showed that very 
distant binary companions may severely affect planetary evo¬ 
lution and influence the orbits of any planet around the other 
component of the system, thus maybe favouring the ingestion 
of m aterial by the host s tar. Fo llowing the results shown in Fig. 
2 of iPinsonneault et~afl ( 2001 1. we roughly estimated that XO- 
2N could have ingested an amount of iron slightly higher than 
5 M® to increase its photospheric iron content by ~0.05 dex, 
given its T c ir ~5300 K. The difference in the iron abundance of 
the two XO-2 stellar components is similar to what found by 
iRamfrez et alJ (1201 ll) for the solar twins 16 Cyg A and 16 Cyg B 
(A[Fe/H] =0.042 ± 0.016 dex). These stars have masses close to 
those of the XO-2 companions and, together with a third com¬ 
panion 16 Cyg C, are members of a hierarchical triple system, 
where the A and C components form a close binary with a pro¬ 
jected separation ~70 AU. The component B is known to host a 
giant planet in a long-period and highly eccentric orbit (P~800 
days and e~O.690. but unlike the case of XO-2 it is the star with 


Mack et al. 2014) did not find relevant differences in elemental 


http://exoplanet.eu 


abundance among the components. All these results imply that 
the presence of giant planets does not necessarily imply differ¬ 
ences in the chemical composition of the host star. 

Our results will produce a more clear picture when several 
iron-peak, a-, s-process, and other odd/even-Z elemental abun¬ 
dances will be analysed for both stars and their behaviour with 
the condensation temperature studied with the aim to investigate 
possible selective accretion of planetary material. One of the ad¬ 
vantages of comparing coeval stars in wide binaries is that any 
difference in their abundance trend could be more related to ac¬ 
cretion of rocky-planetary material rather than to the Galactic 
chemical evolution. This analysis will be the subject of a forth¬ 
coming paper of the GAPS series (Biazzo et al., in prep.), where 
we will investigate if the presence of hot Jupiters could lead the 
host stars to ingest material, which in turn may leave measurable 
chemical imprints in their atmospheric abundances. 

3.3. XO-2 distance and galactic space velocities 

Using our stellar parameters, we derived an estimate of the spec¬ 
troscopic distances of XO-2N and XO-2S by means of the fol¬ 
lowing procedure. We generated Monte-Carlo (MC) normal dis¬ 
tributions for each spectroscopic parameter T e ff, [Fe/H], and 
log#, composed of 10,000 random values and centred on the 
best estimates (Table 0. By keeping the stellar radii fixed to 
the values listed in Table [2] (for XO-2N we used the most ac¬ 
curate estimate), for each MC simulation we first determined the 
stellar bolometric luminosity L * (in solar units) from the Stefan- 
Boltzmann law, and then we derived the absolute bolometric 
magnitude Mbd from the relation Mb 0 i = 4.75 -2.5- log(L»). 
By estimating the appropriate bolometric correction (BC), a 
value for the absolute magnitude in V-band My was then ob¬ 
tained. The BC term was_evaluated using the code provided by 
ICasagrande & VandenBergj ( 2014 ). An additional input is the 
color excess E(B - V ) of the star, which we derived through 
the relation E(B - V) = Ay(s)/3.1, where Ay(s) is the interstel¬ 
lar dust extinction in V-band integrated at the distance s of the 
star (in pc) and measured along the line-of-sight. We derived 
Ay(s) by adopting a simplified model of the local distrib ution 
of the interstellar dust density (’Primmel & Sneraef i200l[) . ex¬ 
pressed by the relation p=po-sech 2 (z/h s ), where z is the height of 
the star above the Galactic plane and li s is the scale-height of the 
dust, for which we adopted the value of 190 pc. The term z is 
related to the distance s and the Galactic latitude of the star b by 
the formula z=s- sin/). From this model we obtained the relation 
Ay(s) = Ay(tot)-sinh(j-- sin/>//z s )/cosh(.s- sin7>//z s ), where Ay(tot) is 
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the interstellar extinction in V-band along the line-of-sight inte¬ 
grated through the Galaxy, and can be estimated from 2-D Galac¬ 
tic maps. For this purpose we used the value A v(tot)=0.16 mag 
derived from the maps of lSchlaflv & Finkbeiner! 0201 ll FT. B y as- 
suming s=150 pc as a prior distance of the stars iBurke et al.l 
1200% . we obtained Av(150pc)~0.06 mag, corresponding to 
E{B - y)=0.019 mag. This is^the value used as input to the code 
of ICasagrande & VandenBergj d2014j) to obtain a first guess of 
the BC in V-band. This in turn was used in the distance modulus 
formula V-(Mb 0 i-BCv)=5- log(s)-5-Ay(s), to obtain a new value 
for the stellar distance s. The new distance was used to repeat the 
procedure iteratively, by determining at each step a new value of 
Ay(s) and BCy(s), and finally another estimate of ,v. When the 
absolute difference between the last and previous calculated val¬ 
ues of s was below 0.1 pc, the iterative process was interrupted 
and the last derived value for s was assumed as the distance of 
the star for the N -th Monte-Carlo simulation. The adopted esti¬ 
mates for the distance of the XO-2 components are the median 
of the distributions of the 10,000 MC values, and the asymmet¬ 
ric error bars defined as the 15.85 ,/j and 68.3'* percentile (see 
Table Q}. Being model-dependent, we do not argue here whether 
the difference of ~2.5 pc between the XO-2S and XO-2S dis¬ 
tances is real. We only note that the two values are compatible 
within the uncertainties and that our best estimate for the dis¬ 
tance of XO-2N locates the star few parsec closer than reported 
bv lBurke et alj d2007l) . 

Distances, together with the stellar equatorial coordinates, 
the proper motions, and the average radial velocities (Table Q] 
and [7}, were used to provide new estimates of the galactic space 
velocities (U,V,W) of the stars with respect to the Local Standard 
of Rest (LSR). We thus assumed (i) a left-handed system of ref¬ 
erence (i.e. the velocity component U is positive toward the di¬ 
rection of the Galactic anti-center), (ii) the transformation matrix 
from equatoria l to galactic coordinates taken from the Hippar- 
cos catalogue dPerrvman & ES All 1 9971) . and (iii) the correction 
for the solar motion in the LSR derived from Coskunoglu et alJ 
(1201 ll) . The results are listed in Table HI [Burke etaL ( 20071) dis- 
cussed the kinematics of the XO-2 system by using data for XO- 
2N only, and concluded that, despite the high galactic space ve¬ 
locity, which indicates a thick disk membership, the super-solar 
metallicity is indeed representative of a thin disk membership. 
This suggests that the high-proper motion of the system can be 
related to an eccentric orbit with a maximum height above the 
galactic plane of ~ 100 pc, then confined to the galactic disk. 

We perf ormed a new analys i s of th e galactic orbits, based on 
the work of lBarbieri & Grattonl (2002), and we confirm the pre¬ 
vious findings. Results for a sample of 1000 orbits show that the 
XO-2 system belongs to the thin disk, with a likely maximum 
height above the galactic plane of 120 pc. This result does not 
change significantly unless the distance of the stars is at least 20 
pc more or less our estimates. Our analysis also indicates a pos¬ 
sible eccentric orbit, with a modal eccentricity equal to 0.42^j® 
(the asymmetric errorbars are defined as the 10'* and 90'* per¬ 
centile), and a minimum distance from the galactic center of 
/?min=3.5±1.5 kpc. If we assume R mm as a rough estimate of the 
galactic region where the two stars formed, this could explain 
their super-solar metallicity, because regions close to the galactic 
center are more metal rich than the solar neighbourhood. More¬ 
over, this could also have modified the orbit of the binary with 
time due to the influence of the more dense stellar environment 
close to the galactic center, with important consequences on the 
evolution of the orbits of the planets. 


available at http://irsa.ipac.caltech.edu/applications/DUST/ 


Table 1 . Basic information about the XO-2 stellar binary system 


Parameter 

XO-2N 

XO-2S 

Ref. 

UCAC4 ID 

702-047113 

702-047114 

(1) 

R.A. (J2000) 

07:48:06.471 

07:48:07.479 

(1) 

DEC (J2000) 

+50:13:32.91 

+50:13:03.25 

(1) 

Near UV (mag) 

18.153±0.047 

17.988+0.044 

(2) 

B (mag) 

12.002+0.015 

11.927+0.015 

(3) 

V (mag) 

11.138+0.026 

11.086+0.025 

(3) 

Rc (mag) 

10.669+0.020 

10.631+0.019 

(3) 

L (mag) 

10.243+0.012 

10.216+0.013 

(3) 

J (mag) 

9.744+0.022 

9.742+0.022 

(4) 

H (mag) 

9.340+0.026 

9.371+0.027 

(4) 

K, (mag) 

9.308+0.021 

9.272+0.021 

(4) 

Distance (pc) 

145.511;® 

148.0111 

(5) 

Ba (mas yr 1 ) 

-31.9+3.0 

-29.9+3.6 

(1) 

Bs (mas yr 1 ) 

-156.2+3.3 

-156.2+3.6 

(1) 

U (km s _1 ) 

70 39 +0 - 79 

71.051” 

(5) 

V (km s _1 ) 

-76 0 +2 ° 

/u - 2.0 

-77.651J5 

(5) 

W (km s _1 ) 

—9 Q 1+0.67 

1 -0.66 

-3-46tg 

(5) 

Relative magnitudes (XO-2N - XO-2S) 

AB (mag) 

0.075+0.002 

(6) 

AV (mag) 

0.057+0.002 

(6) 

AR ( . (mag) 

0.040+0.002 

(6) 


0.041+0.004 

(7) 

Ar' (mag) 

0.038+0.001 

(8) 

AI C (mag) 

0.031+0.003 

(6) 


0.030+0.006 

(9) 


References. (1) UCAC4 catalogue dZacharias et all 1 201 3l): (2) 
iBianchi et all d2011l) : (3) APA SS all-sky survey iHenden & Munaril 
d2QI4ll : (41 2MASS catalogue dSkrutskie et al.ll2006h : (51 This work: 
(6) This wor k (Serra La Nav e Observatory); (7) This work (TASTE 
project); (8) dKundurthv et al.l[20 1 3l) . Ai,~626 nm, 4 epochs ; (9) This 
work (APACHE project). 


4. Analysis of the light curves 

4. 1. Out-of-transit photometry 

Figure [4] shows the relative photometry of the two components 
of the XO-2 system obtained from different facilities. Measure¬ 
ments were collected in the BVRI standard photometric pass- 
bands and the points corresponding to the transits of the planet 
XO-2Nb were not included in the calculation of the intra-night 
average values. The uncertainties associated to the measure¬ 
ments are the RMS of the points collected in each single night 
of observation, and when only one measurement per night was 
available from SLN, we adopted the dispersion of all the data 
as the error. The average values for the magnitude differences in 
each band are indicated in Table |T] We note that the indepen¬ 
dent measurements in R (TASTE and SLN) and I (APACHE and 
SLN) bands agree well each other. Our derived relative magni¬ 
tudes have been used in the differential abundance analysis dis¬ 
cussed in Sect l3.2l 
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Table 2. Stellar parameters for the XO-2 components derived from the 
analysis of the HARPS-N spectra and from stellar evolutionary tracks. 
The adopted values represent the weighted mean of individual measure¬ 
ments obtained with different methods (except for the projected velocity 
Vsin/*, for which we adopt the measurement derived from the Rossiter- 
McLaughlin effect), and their associated uncertainties are, for a more 
conservative estimate, the average of the uncertainties of the individual 
values. 


Parameter 

XO-2N 

XO-2S 

Note 

T e// [K] 

5332+57 

5395+54 


log g [cgs] 

4.44+0.08 

4.43+0.08 


[Fe/H] [dex] 

0.43+0.05 

0.39+0.05 


Microturb. £ [km s' 1 ] 

0.88+0.11 

0.90+0.10 


Vsin/* [km s' 1 ] 

1.07+0.09 

1.5+0.3 


Mass [M e ] 

0.97+0.05 

0.98+0.05 

(1) 


0.96+0.05 

- 

(2) 

Radius [Rq] 

1 01 +01 
1 - Ui -0.07 

1 02 +009 

1-U -0.06 

(1) 


0 998 +0 033 

U.7!7O_0 032 

- 

(2) 

Age [Gyr] 

7 0+2.3 
'• -3.0 

7 1 +2 - 5 

'-2.9 

(1) 


7-8Tl;l 

- 

(2) 

Luminosity [L 0 ] 

0.70+0.04 

0.79+0.14 

for XO-2N: (2) 


for XO-2S: (1) 


Notes. (1) Matching the T eff , [Fe/H], and \ogg to the Yonsei-Yale evo¬ 
lutionary tracks. (2) Matching the T e fr, [Fe/H], and stellar density to the 
Yonsei-Yale evolutionary tracks. 


Table 3. Stellar atmospheric parameters and their differences for XO- 
2N and XO-2S. ee discussion in Sect l3.2l 


Parameter 

XO-2N 

XO-2S 

An-s 

Analysis based on equivalent widths 

Teff [K] 

5320+50 

5330+50 

-10 + 71 

log g [cgs] 

4.46+0.08 

4.41+0.12 

+0.05 + 0.14 

£ [km s' 1 ] 

0.89+0.14 

0.93+0.04 

-0.04 + 0.15 

[Fe/H] fl [dex] 

0.37+0.07 

0.32+0.08 

+0.05 + 0.11 


Line-by-line analysis 


Teff [K] 

5290+18 

5322+25 

-32 + 31 

log g [cgs] 

4.43+0.10 

4.41+0.10 

+0.02 + 0.14 

£ [km s' 1 ] 

0.86+0.06 

0.93+0.05 

-0.07 ± 0.08 

[Fe/H] fl [dex] 

0.37+0.07 

0.32+0.08 

+0.05 + 0.11 


Differential analysis 


Teff (K] 

5290+18 

5325+37 

-35 + 8 

log g [cgs] 

4.43+0.10 

4.420+0.094 

+0.010 + 0.020 

£ [km s' 1 ] 

0.86+0.06 

0.93+0.03 

-0.07 ± 0.07 

[Fe/H]" [dex] 

0.37+0.07 

0.32+0.08 

+0.054 + 0.013 


Notes. (u) The iron abundance [Fe/H] refers to the abundance of Fe i. 


On the night between October 31, and November 1 2014, 
we took for 30 minutes a series of 60 s exposures of the XO- 
2 field with the Asiago Faint Object Spectrograph and Camera 
(AFOSC) at the 182cm Asiago Astrophysical Observatory tele¬ 
scope, while the field was close to the zenith. By placing both 
components in a 8" slit and using a grism covering the spectral 
range 3900-8000 A, we obtained a final simultaneous differential 
spectrum by dividing each other the combined individual spectra 
in each wavelength resolution element. The differential spectrum 
is shown in Fig. 0 These observations were characterised by a 



X (eV) 


Fig. 3. Iron abundance difference between XO-2N and XO-2S versus x- 
Dashed line shows the mean A[Fei/H], while dotted lines represent the 
standard deviation ±lcr from the mean, where <+=0.022 dex. 


very good instrumental stability and turned out to be useful to 
strengthen the results obtained with broad-band measurements 
of the differential magnitudes. In fact, these are in good agree¬ 
ment with the AFOSC measurements, as shown in Fig. [5] Due to 
the high scatter of the data in the /-band caused by fringing, we 
did not include in the plot the measurements with wavelengths 
greater than ~770 nm. 

4.1.1. Rotational periods 

Thanks to the timespan of the APACHE /-band observations and 
~800 useful data points for each component, we could analyse 
the photometric time series to search for periodic modulations 
ascribable to the rotation of the stars and produced by brightness 
inhomogeneities distributed on the stellar photospheres (which 
are assumed to change on time scales longer than the timespan 
of the observations). Taking into account that the XO-2 system 
is quite old (see Table [2} and the two stars are magnetically quiet 
(see discussion in Sectj6]», they are expected to have rotation pe¬ 
riods of some tens of days and low amplitudes in the light curve 
variations^ We_used the Generalized Lomb-Scargle algorithm 
(GLS: IZechmeister & Kiirsteil12009 ) to search for periodicities 
in our photometric time series, after binning in a single point the 
data collected during each night. The data are available on-line 
at the CDS, each containing the timestamps, differential magni¬ 
tudes and corresponding uncertainties. Results are shown in Fig. 
[6]and[7]for XO-2N and XO-2S. The major peaks at low frequen¬ 
cies in the periodograms are due to the limited time baseline of 
our observations. For the North component we find a peak in the 
periodogramcorresponding to a period P=41.6±l.l days, with a 
semi-amplitude of the folded light curve of 0.0027+0.0004 mag. 
Because there is no evidence of harmonics in the periodogram, 
the measurements can be arranged almost perfectly with a sine 
function. To estimate the significance of the best period, we per¬ 
formed a bootstrap analysis (with replacement) by using 10,000 
fake datasets derived from the original photometry, and in six 
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cases we found a power greater than that associated to the best 
period in the real dataset. This corresponds to a False Alarm 
Probability (FAP) equal to 0.06%, which makes the result for 
XO-2N rather robust. 

Interestingly, to our knowledge other two systems with an 
old star have measured rotation periods thanks to Kepler data. 
The A and B components of the 16 Cygni system, which is 
estimated to be 6.87+2 Gyr old, have rotation periods 23.8' ['jj 
and 23.2^2 days respectively ( iDavies et al.l 2015 ). Another 
example is that of Kepler-77, a star-planet system 7+2 Gyr 
old. The detected stellar rotation period is 36+6 days, with 
sizeable changes in the shape of the light curve modulation 
(iGandolfi et alJl2013 ). This case and especially that of 16 Cyg 
B, with uncertainties of several days, demonstrate that for old 
late-type stars the determination of the true rotation period can 
be particularly difficult. The intrinsic evolution of their surface 
active regions takes place on a time scale shorter than the rotation 
period, and therefore starspots become bad tracers of stellar ro¬ 
tation. In the Sun, depending on the phase of the solar cycle, the 
rotation period measured by fitting a sinusoid to the total solar 
irradiance - a good proxy for the Sun-as-a-star optical modula¬ 
tion - varies between 22.0 and 31.9 days (the true synodic period 
being 27.27 days). The variations are especially large during the 
maximum of the solar cycle because several large active regions 
are simultaneously present on the solar surface and evolve on 
time scales of 1-2 weeks, i.e. significantly shorter than the true 
rotation period (e.g., lLanza et all2003}) . 

Even if the rotation period of XO-2N appears well con¬ 
strained, however the complications expected for old late-type 
stars could actually have prevented us to detect an accurate value. 
This is indeed the case for the companion star XO-2S, for which 
a less convincing detection resulted in our data. The GLS peri- 
odogram returned the highest peak signal at / =0.0384 days -1 , 
corresponding to P=26.0±0.6 days, with a semi-amplitude of 
0.0024+0.0005 mag of the folded light curve. Also for XO-2S 
there is no evidence of harmonics. The bootstrap analysis in this 
case led to a LAP=2.6%, which casts some doubts on the reality 
of the signal, even if the period is still compatible with the mea¬ 
sured projected velocity Vsin/*. The similarity between the ob¬ 
served periodogram and the recentered spectral window (bottom 
panel in the upper part of Fig. |7]i strengthens our confidence on 
the stellar origin of the periodicity we found. However, it must 
be noticed a second peak at 0.029 days -1 (P=34.5 days) with a 
slightly less power. This has some relevance in the discussion, 
since a rotational period of 34.5 da ys mat ches the gyrochronol- 
ogy relations (e.g. see Fig. 13 in lBarnesl2007l) better than the one 
of 26 days, and it is also closer to the estimates based on mea¬ 
surements of the activity index R' HK (see Table [8]). Moreover, in 
analogy with the case of the Sun-as-a-star, if we assume P=34.5 
days as the best detection, a true period which is a bit higher than 
that is nonetheless reliable. Therefore more data are necessary, 
possibly covering several activity cycles, to reach a more robust 
conclusion concerning the true rotational period of XO-2S. 


Through the combined analysis of 14 new transit light 
curves, we present here refined parameters for the XO-2N sys¬ 
tem, together with those in Table [2] The individual transit light 
curves are shown in Fig. [8j and the full dataset combined in 
a single, high-quality transit light curve is shown in Fig. [9] 
The information about each dataset is summarized in Table [4] 
The raw scientific frames were first corrected for bias, dark and 
flat-fi eld using standard softw are tools. The STARSKY pipeline 
(see iNascimbeni et al.l 1201 3l and references therein) was em¬ 
ployed to perform differential photometry between XO-2N and 
a set of optimally-weighted reference stars. The whole set of 
fourteen light curves was then simultaneously analysed with 
the JKTEBOfQ code (version 34, ISouthworth et alJ 12004). The 
adopted transit model parametrizes the stellar limb darkening 
(LD) with a quadratic law, i.e. with the two parameters u\ and 
ii 2 - Assuming a circular orbit, seven parameters were left free 
to vary in the fitting procedure: the period P and reference time 
Tq of the underlying transit ephemeris, the ratio k-Rp/R^ and 
sum (Rp+R^/a of the fractional radii (a being the semi-major 
axis of the system), the orbital inclination i, and the LD param¬ 
eters u i and i< 2 - Leaving both u\ and in free to vary allows us to 
avoid to put any a priori assumption on the stellar atmospheric 
parameters. Our estimates are in good agreement with those de- 
rivedby interpolating the theoretical tables of lClaret & Bloemenl 
1 2011 ) using our derived stellar atmospheric parameters (0.48 
and 0.21 for u\ and ui). The transit model was fitted to our data 
with a non-linear least-squares technique; the uncertainties over 
the best-fit parameters are then recovered with a bootstrap al¬ 
gorithm, described by ISouthworthl d2008l) . Other quantities of 
interest, such as the impact parameter b or the transit duration 
7 ) 4 , were derived through their analytical expressions from the 
fitted parameters. The final best-fit values, along with their esti¬ 
mated errors, are summarized in Table 0 We give the estimates 
as the median of the bootstrapped distributions, with the asym¬ 
metric errorbars defined as the 15.87 ,/! and 84.1 3 th percentile. 
Our results are in good agreement with those in the literature, in 
particular with the results of ICrouzet et al.l 12012)) based on ob¬ 
servations made with the Hubble Space Telescope. We were able 
to provide an updated ephemeris (with an error on P of about 
0.025 s) that we also used as a prior for the radial velocity fit 
(Sect]5j. Our determination of the stellar density p/p G , thanks to 
u i and ui being left free to vary, is completely independent from 
stellar models, and was used as input to interpolate within the 
Y-Y evolutionary tracks and derive an accurate estimate for the 
stellar mass, radius and age (Sect0. 

5. Radial velocities of XO-2N 

We observed the star XO-2N as part of a program aimed at 
searching for additional outer planets in the system, and we in¬ 
vestigated the Rossiter-McLaughlin (RML) effect by acquiring 
high-resolution spectra with HARPS-N during one transit event 
of the planet XO-2Nb. 


4.2. Transit photometry of XO-2N 

Lollowing the discovery of the hot Jupiter XO-2Nb dBurke et alJ 
12007I) . new transit light curves of XO-2N were collected to im¬ 
prove the stellar a nd planetary parameters and search for transit 


Crouzet et al.l 12012 

). and the planet atmosphere was analysed 

with space-based 

Machalek et al. 2009; 

Crouzet et al. 2012) 

and ground-based observations (Sing et al.1120 1 2012b- 


5.1. Analysis of the radial velocity time series 

To refine the orbital parameters of XO-2Nb, our new HARPS- 
N data were analysed along with the RVs obtained with the 
spectrographs HDS (High Dispers ion Spectrograph) at the Sub¬ 
aru telescope dNarita et alJ 1201 lT) . and HIRES (High Resolu¬ 
tion E chelle Spectrograph) at the Keck teles cone dKnutson et al.l 
i2014l) . The original spectroscopic data from lBurke et al.l ( 2007!) 


http://www.astro.keele.ac.uk/jkt/codes/jktebop.html 
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Table 4. Log of the transit observations of XO-2Nb analysed in this work (14 distinct light curves). All of them were obtained using a 
Cousin/Bessell R filter. 


Date 

Observatory 

Telescope/instrument 

N 

cr 

(mmag) 

cr i20j 
(mmag)) 

2011/01/03 

Asiago 

182cm/AFOSC 

2620 

1.37 

0.37 

2011 /02/01 

Asiago 

182cm/AFOSC 

2483 

2.19 

0.42 

2011/02/06 

Asiago 

67/91cm Schmidt 

1029 

3.98 

1.35 

2011/02/06 

Asiago 

182cm/AFOSC 

2615 

2.16 

0.51 

2011/04/02 

Asiago 

182cm/AFOSC 

1767 

2.41 

0.52 

2011/11/26 

Asiago 

182cm/AFOSC 

2999 

2.28 

0.50 

2011/12/17 

Asiago 

182cm/AFOSC 

2343 

1.59 

0.36 

2012 /01/12 

Teide 

IAC80/CAMELOT 

899 

2.19 

0.86 

2012/02/23 

Teide 

IAC80/CAMELOT 

955 

1.77 

0.55 

2012/03/15 

Asiago 

182cm/AFOSC 

2853 

1.63 

0.47 

2012/12/05 

Teide 

IAC80/TCP 

712 

1.85 

0.66 

2012 /12/12 

Asiago 

182cm/AFOSC 

2736 

1.41 

0.37 

2013/02/05 

Asiago 

182cm/AFOSC 

1277 

2.38 

0.66 

2014/03/07 

Asiago 

182cm/AFOSC 

2114 

2.19 

0.59 


Notes. The columns list: the evening date, the observatory and instrument employed, the number of data points N, the photometric RMS of the 
residuals from the best-fit model (both for the original cadence, cr, and after binning over 120-s intervals, cr 120 ). 

Table 6. Orbital parameters for the XO-2N and XO-2S systems. Errors and upper limits_refer to ltr uncertainties. The results for XO-2Nb are 
discussed in Sect. EH while those for the XO-2S planets are taken from lPesidera et a~l. 12014 ). We also include our estimates of few parameters 
for a second companion to XO-2N, namely XO-2Nc, despite its existence is a matter of debate, as discussed in Sect. l5.lTl 


Parameter 

XO-2Nb 

XO-2NT 

XO-2Sb 

XO-2Sc 

Note 

P [days] 

2.61585922 ± 2.8 • 10~ 7 

>17 years 

18.157±0.034 

120.80±0.34 

fitted 

K [m s~' ] 

90.17 ±0.82 

> 20 

20.64±0.85 

57.68±0.69 

fitted 

K ■ P~ 2 [ms _1 day~ 2 ] 

- 

-5.13 ± 0.54 • 10~ 7 

- 

- 

derived 

-\/esincu 

- 

- 

-0-314!^ 

U. Joo_o 012 

fitted 

yfecoso) 

- 

- 

0.282^ 

-0.0381“ 33 

fitted 

e COS Cl) 

-0 0009 + ° 0034 

U.UUU7_ 00022 

- 

- 

- 

fitted 

e sintd 

-0 0004 + ° 0040 

- 

- 

- 

fitted 

e 

< 0.006 

- 

0.180 ±0.035 

0 1 +0.0094 

u.i JZ,o_ 0009g 

derived 

0 ) [deg] 

unconstrained 

- 

311.9 ± 9.5 

264.5 ± 4.9 

derived 

T c [BJD tdb -2,450,000] 

5565.546480 ± 5.4 ■ 10- 5 

- 

6413 11 +0 - 82 

uh-id.ii_ 086 

6408.1/| p 

fitted 

ThdsDus -1 ] 

-17.3 ± 2.0 

- 

- 

- 

fitted 

"Xhires [ms '] 

-13.5 ±2.9 

- 

- 

- 

fitted 

Lharps-n [ rn s '] 

46932.9 ± 1.7 

- 

46543±1 

fitted 

;zffer H Ds[ms~ l ] 

< 1.2 

- 

- 

- 

fitted 

yzfter H iREs[ms -1 ] 

7 0 +2 ' 9 
' -2.0 

- 

- 

- 

fitted 

/;tter H ARps-N[ m s -1 ] 

2.6 ± 0.7 

- 

1.80 ±0.43 

fitted 

y [ m s~*day _1 ] 

0.0017 ± 0.0018 

- 

0.0531±0.0087 

fitted 

y [ms -1 day -2 ] 

-2.03 ± 0.21 ■ 10~ 5 

- 

- 

- 

fitted 

M sin i [My„ p ] 

0.597 ± 0.021“ 

> 1.8 

0.259 ± 0.014 

1.370 ±0.053 

derived 

Radius R [Rj up ] 

1.019 ±0.031 

- 

- 

- 

derived 

Density p [g cm -3 ] 

0 70+ 0072 
' u -0.064 

- 

- 

- 

derived 

Surface gravity log g [cgs] 

3.15 ±0.03 

- 

- 

- 

derived 

a [AU] 

0.03673±0.00064 

- 

0.1344 ±0.0025 

0.4756 ± 0.0087 

derived 


Notes. (o) For XO-2Nb this corresponds to the real mass of the planet. It was obtained with a Monte-Carlo analysis using the median values of the 
orbital period and inclination angle of the orbital plane listed in Table El the stellar mass in Table [2] (second value), and the amplitude K of the 
radial velocity included in this Table. 
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Table 7. Radial velocities (RV), bisector velocity span (BIS) and activity index logfR^) of XO-2N, as measured with the HARPS-N spectrograph. 
The BIS reported here are those derived by the data reduction pipeline of HARPS-N, with uncertainties assumed to be twice those on the radial 
velocities. 


Time 

(BJD UTC - 2,450,000) 

RV 

(km s' 1 ) 

RV error 
(km s' 1 ) 

BIS 

(km s' 1 ) 

iogl^HK) 

error log(R^ K ) 

6252.773594 

47.0169 

0.0047 

0.0209 

-5.172 

0.157 

6266.500236 

46.9521 

0.0018 

0.0204 

-4.871 

0.022 

6266.511139 

46.9500 

0.0016 

0.0213 

-4.870 

0.019 

6266.522077 

46.9472 

0.0019 

0.0187 

-4.905 

0.026 

6266.533015 

46.9416 

0.0018 

0.0164 

-4.907 

0.025 

6266.543953 

46.9406 

0.0021 

0.0210 

-4.927 

0.033 

6266.554891 

46.9490 

0.0028 

0.0303 

-4.874 

0.043 

6266.565829 

46.9422 

0.0024 

0.0245 

-4.891 

0.035 

6266.576768 

46.9413 

0.0028 

0.0249 

-4.986 

0.054 

6266.587706 

46.9388 

0.0023 

0.0150 

-4.889 

0.034 

6266.598644 

46.9305 

0.0023 

0.0229 

-4.855 

0.031 

6266.609593 

46.9217 

0.0022 

0.0274 

-4.881 

0.030 

6266.620520 

46.9181 

0.0018 

0.0195 

-4.889 

0.023 

6266.631458 

46.9153 

0.0017 

0.0205 

-4.910 

0.023 

6266.642396 

46.9162 

0.0016 

0.0178 

-4.901 

0.019 

6266.653334 

46.9196 

0.0015 

0.0209 

-4.872 

0.017 

6266.664284 

46.9179 

0.0016 

0.0217 

-4.901 

0.020 

6266.675210 

46.9135 

0.0015 

0.0191 

-4.890 

0.017 

6266.686148 

46.9138 

0.0015 

0.0152 

-4.877 

0.017 

6297.623295 

46.9958 

0.0015 

0.0154 

-4.928 

0.019 

6298.643938 

46.8409 

0.0022 

0.0097 

-4.921 

0.034 

6299.762797 

47.0128 

0.0015 

0.0099 

-4.907 

0.020 

6305.631219 

46.9708 

0.0020 

0.0125 

-4.938 

0.031 

6322.563804 

46.8777 

0.0015 

0.0129 

-4.898 

0.019 

6323.686453 

47.0147 

0.0033 

0.0209 

-4.947 

0.072 

6380.349341 

46.9211 

0.0033 

0.0106 

-5.269 

0.155 

6400.444868 

46.8469 

0.0100 

0.0573 

-4.861 

0.214 

6406.361875 

46.8938 

0.0036 

0.0113 

-4.908 

0.060 

6420.392578 

47.0089 

0.0031 

0.0066 

-5.229 

0.114 

6437.372927 

46.8412 

0.0014 

0.0140 

-4.964 

0.018 

6548.731988 

47.0043 

0.0024 

0.0059 

-4.949 

0.044 

6586.684893 

46.8605 

0.0024 

0.0145 

-4.989 

0.044 

6602.751791 

46.9259 

0.0025 

0.0094 

-5.042 

0.055 

6616.709751 

47.0041 

0.0027 

0.0203 

-4.874 

0.042 

6617.722042 

46.8368 

0.0017 

0.0153 

-4.929 

0.023 

6631.704304 

46.9647 

0.0020 

0.0051 

-5.035 

0.041 

6657.453439 

46.8928 

0.0056 

0.0255 

-5.086 

0.183 

6700.404020 

47.0066 

0.0026 

0.0147 

-4.983 

0.051 

6701.349611 

46.8410 

0.0019 

0.0164 

-4.963 

0.031 

6783.385415 

46.9513 

0.0014 

0.0077 

-4.996 

0.024 

6784.360848 

46.9543 

0.0016 

0.0082 

-4.980 

0.028 

6933.733982 

46.9031 

0.0023 

0.0052 

-5.009 

0.050 

6935.754527 

47.0090 

0.0017 

0.0107 

-4.996 

0.030 


were not included because they have considerably larger error 
bars (> 20 ms -1 ) than the other datasets and thus do not yield 
any significant improvement of the orbital solution. We first fit¬ 
ted separately the HDS, HIRES, and HARPS-N RV time series 
with a Keplerian orbit and a slope, after removing the in-transit 
points which are affected by the Rossiter-McLaughlin effect. 
We then performed a Bayesian analysis with our Differential 
Evo lution Markov Chain Monte Carlo (DE- MC MC) tool (see, 
e.gJDesidera et al.l2014l) . lNarita et all (1201 ll) and lKnutson et al.l 
(1201411 reported on a positive acceleration revealed in their RV 
datasets, equal to 0.0206 ± 0.0016 and 0.0126^^^ m s'May' 1 
respectively, possibly due to a second massive companion in an 
outer orbit, namely XO-2N c. Our observations and analysis con¬ 
firm the existence of this trend, but we observed a turn-over in 
the HARPS-N data and measured a negative acceleration equal 
to -0.0170!°““ ms _1 day _I . 


We then modelled the RV data of the three datasets with a 
Keplerian orbit and a polyn omial model following the formal¬ 
ism of iKinning et al.l «201 li) and choosing as f plvol in their Eq. 
(1) the mean of the epochs of the RV observations. In total, 
our RV model has thirteen free parameters: the transit epoch 
T c , the orbital period P, the radial-velocity semi-amplitude K, 
e cos ix) and e sin at of XO-2N b, where e and at are the eccentric¬ 
ity and argument of periastron; three systemic velocities (thds, 
Thires, Tharps-n) and RV jitter terms (jitter HD s, y/ffer H iRES- 
jitteruARPS-N), corresponding to the three different datasets; the 
linear f and quadratic y trends. Gaussian priors were imposed 
on i) T c and P from the new ephemeris we derived in Sect]?] 
and ii) the centre times of the secondary eclipse observed by 
iMachalek et aD (i2009l) at 3.6, 4.5, 5.8, and 8.0 pm with Spitzer, 
because these provide strong constraints on the orbital eccen¬ 
tricity; for this purpose, we used Eq. (18) in IJordan & Bakosl 


Article number, page 12 of [26] 























M. Damasso et al.: GAPS V: Global analysis of the XO-2 system 



Fig. 5. Upper panel. Differential spectro-photometry of the XO-2 com¬ 
ponents obtained with the AFOSC camera at the 182cm Asiago tele¬ 
scope. Red dots represent the mean values of measurements in bin of 
200 A, while blue squares indicate the broad-band differential magni¬ 
tudes presented in Table |T] Lower panel. Bandpasses of the filters used 
for the broad-band measurements. 

Table 5. Transit parameters derived from the combined analysis of the 
14 transit light curves of XO-2Nb plotted in Fig. [ 8 ] 


Parameter 

Value 

P(d) 

2.61585922+0.00000028 

(R p + R*)/fl 

0.1393+0.0018 

R P /R+ 

0 10490 + ° 00059 

U.1UH-!7U_ 000063 

(R P /R*) 2 

0.01100+0.00013 

i (deg) 

87.96 

a/R* 

7 Q9Q+0-099 

1 -0.093 

b 

o-2S7-;;;;:i 

T 14 (min) 

162.15+0.36 

T 14 (d) 

0.11260+0.00025 

T c (BJD TDB-2,450,000) 

5565.546480+0.000054 

U\ 

0 474 + 0.030 
-0.028 

U2 

0 171+ 0067 
' 1 -0.070 

P/Pe 

0 960 +0037 

U,7U -0.034 


Notes. The listed parameters are: orbital period P, sum and ratio of 
the fractional radii (Rp+R*)/a, (Rp/R*), orbital inclination i, semi-major 
axis a scaled by stellar radius R*, impact parameter b, total duration of 
the transit from the first to the fourth contact 7 ) 4 , reference time of the 
ephemeris T 0 (time standard is Barycentric Julian Day calculated from 
UTC times: BJDutc), the linear and quadratic LD parameters ui , 112 , the 
average stellar density p/p B . 

(l2008l). We adopted Jeffrey’s priors for the three jitter terms (see, 
e.g.. lGregorvl2005l) and non-informative priors for the RV semi¬ 
amplitude and the three systemic velocities. 

After the earliest analysis, we noticed that the first HIRES 
point at 2454372.14 BJDutc was 3cr below our best fit hence 
we discarded it and repeated the analysis, although its inclu¬ 
sion has a very minor influence on the final solution. The pos¬ 
terior distributions of the fitted and derived parameters were ob¬ 
tained with the same procedure as described in iDesidera et ali 


APACHE photometry periodogram of the star X0-2N 



APACHE photometry of the star X0-2N 



Fig. 6. Upper panel. The GLS periodogram of the XO-2N light curve 
collected by the APACHE survey and the corresponding window func¬ 
tion, with the highest peak translated to the frequency of the best signif¬ 
icant peak found in the periodogram. Levels of different FAP are also 
indicated. Lower panel. The APACHE light curve of XO-2N folded at 
the best period found by GLS with superposed the best sinusoidal fit of 
semi-amplitude A=0.0027 mag (solid red line). 

( 2014 ). Their medians and 68.3% confidence intervals are re¬ 
ported in Table [6] along with the system parameters of XO-2S 
derived bv iDesidera et al.l (i2014 L The Keplerian best fit to the 
three datasets, after removing the quadratic trend, is shown in 
Fig. [TO] Our new RV semi-amplitude K = 89.78 + 0.81 m s" 1 is 
slightly lower than the value found bv lKnutson et ali (2014), that 
is K = 93.9 ± 2.2 m s _1 . Its precision is improved by a factor of 
three, especially thanks to the HARPS-N data. The orbital eccen¬ 
tricity is consistent with zero, e < 0.006 at lcr, which justifies the 
use of a circular model for the transit modelling (Sect l4.2l >. The 
combined analysis of the complete RV dataset and of the new 
transit light curves allows us to refine the physical parameters of 
the planet XO-2Nb, and we find that the best values for the ra¬ 
dius, mass and density are 1.019 + 0.03lRj up , 0.597±0.021Mj up , 
and 0 . 70 ™ g cm -3 respectively. 

5.1.1. RV acceleration in XO-2N: evidence for an outer 
companion or solar-like stellar activity cycle? 

Under the hypothesis that an outer massive companion exists and 
is responsible for the observed acceleration, the change of sign 
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APACHE photometry periodogram of the star XO — 2S 



0.0 0.1 0.2 0.3 


Frequency (1 /d) 


APACHE photometry of the star X0-2S 



Fig. 7. Upper panel. The GLS periodogram of the XO-2S light curve 
collected by the APACHE survey and the corresponding window func¬ 
tion, with the highest peak translated to the frequency of the best signif¬ 
icant peak found in the periodogram. Levels of different FAP are also 
indicated. Lower panel. The APACHE light curve of XO-2S folded at 
the best period found by GLS with superposed the best sinusoidal fit of 
semi-amplitude A=0.0024 mag (solid red line). 


derived from HARPS-N data with respect to the previous mea¬ 
surements could be likely caused by the curvature of the orbit 
of XO-2N c. Figure QT| displays the quadratic trend in the RV 
residuals, after removing the orbital signal of XO-2Nb and con¬ 
sidering all the available datasets. By assuming a circular orbit 
for the outer companion XO-2N c, the derived quadratic trend 
allows us_to estimate its orbital period and minimum mass (see 
Kipp ing et al.ll201 H) : P > 17 yrs and Msinf >1.8 Mj up . Direct 
imaging observations collected with the AstraLux instrument did 
not reveal a companion to XO-2N with A;' < 8 magnitudes and 
a separation < 2" (iDaemgen et al .1 [20091) . which corresponds to 
~300 AU at the star location. A rough estimate for an upper limit 
of the XO-2N c mass can be made using the information on A/'. 
By fixing the age and the metallicity to the values of XO-2N, and 
using the CME(_I stellar isochrones, we find that the unseen com¬ 
panion should have a mass M<0.16 M®, which roughly corre¬ 
sponds to a main-sequence star of spectral type later than M5V- 
M6V. 

If not due to an outer companion, the quadratic trend could 
reflect a solar-like activity cycle, which leaves an imprint in the 



0.02 0 0.02 
orbital phase 


Fig. 9. Photometric data presented in Fig.[8]and combined together in a 
single transit light curve to appreciate the quality of the whole dataset. 
Data points are the weighted mean of the measurements in bins of 120 
seconds, and the best-fit transit model is over plotted as a solid line, with 
the O-C residuals shown in the lower panel. 


RV measurements over a timespan of several years. This is con¬ 
sistent with the observation that older stars, with longer rota¬ 
tion periods, have l onger activity c ycles (e.g., Noyes_etaL|[l984; 
IOssendriiverfll997l) . In their work. [Knutson et al.l (1201 4l) do not 
mention for XO-2N the existence of a correlation between their 
RV residuals and some indicators of the stellar activity, as R' HK 
or .S'hk, which would call the planetary origin of the acceleration 
into question. The same situation does not occur in our dataset. 
The time series of the R uy activity index, calculated according to 
the recipe described in lLovis et al.l (1201 111 , shows the existence 
of a negative trend (see Fig. S3. We modelled the data with a 
linear fit, obtaining a negative slope with a lcr uncertainty of 
~ 14%. The Spearman’s rank correlation coefficient of the data 
is quite high (p=-0.79) and with a low FAP (1.3TO~ 6 0 thus 
strengthening the real existence of a correlation. The direction 
of the trend being similar for the RV residuals and R' HK , a sig¬ 
nificant positive correlation is observed between them (see Fig. 
US. characterised by a Spearman’s rank correlation coefficient 
p=+0.68 (FAP = 1-10~ 4 ) forN=26 measurements. 

Due to the importance of this result, the significance of 
the correlation was further investigated by applying two differ¬ 
ent tests. Through the first method, described in IFigueira et al.l 
(12014l) (and references therein) to which we refer to as the null 
hypothesis approach, we evaluated the probability of having a 
Spearman’s rank correlation coefficient larger or equal than that 
of the real data by testing the null hypothesis, i.e. by assum¬ 
ing that the data pairs are actually uncorrelated. This method, 
which does not take into account the data uncertainties, employs 
a Fisher-Yates data shuffling to generate fake data pairs which 
are uncorrelated. We generated 10,000 such random datasets 
and then calculated the z-score for the correlation coefficient 
of the original dataset relative to the distribution of the simu- 


9 http://stev.oapd.inaf.it/cgi-bin/cmd 


10 Values obtained with the R_C0RRELATE function in IDL. 
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lated Spearman’s rank correlation coefficients. The correspond¬ 
ing probability (the p-value) is calculated as the one-sided like¬ 
lihood of having such a z-score from the observed Gaussian dis¬ 
tribution. The second method, which instead takes into account 
the data uncertainties, is described in lCurra 3 (120141 ) (referred to 
as composite method) and it employs a double Monte-Carlo re¬ 
sampling of the dataset pairs. First, fake data pairs are generated 
through a bootstrap (with re-sampling) method, with new data 
pairs randomly chosen from the original ones (i.e. each new pair 
comes from the real dataset and it may appear more than once 
or not at all in a simulated dataset). Then, the second step con¬ 
sists in randomly shuffling each measurement within its uncer¬ 
tainty following a normal distribution of width 1 and centred on 
0. These two steps are repeated 10,000 times, resulting in a dis¬ 
tribution of Spearman’s rank correlation coefficients for which 
the mean and the standard deviation are calculated. Moreover, 
the z-score is calculated for each simulated dataset through the 
Fisher transformation of the correlation coefficient, as discussed 
in ICurranl d2014l) . resulting in a z-score distribution for which 
the mean and the standard deviation are provided. The z-score 
defines the level of significance of a correlation, being the num¬ 
ber of cr from a null correlation for the original dataset. 

The null hypothesis test applied to the RV residuals and R' HK 
dataset for XO-2N resulted in z-score = -3.41, corresponding 
to p-value = 0.032%, while the composite method resulted in 
Pspearman-0-54-±0A5 and z-score = 1.267+0.451. Both results 
suggest that the null-hypothesis can be rejected for our dataset. 

A linear regression applied to the same dataset results in a 
slope of 30.2 m s 1 pci' unit of with a lcr uncertainty of 
~ 10%. Comparing this slope value with that_predicted by the 
empirical model described by Eq. (13) in lLovis et alJ (i201 ll) for 
an inactive star with the same 7' c (r and [Fe/H] of XO-2N, i.e. 
32.9+0.2 ms~'per unit of R' HK , we find that our estimate is in 
agreement with that expected for RV variations produced by a 
long-term stellar activity cycle. 

The positive correlation found between the RV residuals and 
/(J |K indicates that the variations seen in the RV residuals are 
likely related to the intrinsic behaviour of the star. In fact, dur¬ 
ing a phase of high stellar activity a net redshift should appear in 
the RV measurements as a consequence of the suppressed con¬ 
vective motion in the photospheric activ e regions, which gives a 
blueshift contribution (see, for instance, Dumus aue et akl 20121 
and references therein). This means that a negative trend in the 
RVs is expected when a star is observed during a declining phase 
of its stellar cycle, with a similar trend appearing in the R' HK in¬ 
dex. Despite in our data there is evidence in favour of the stellar 
activity cycle as the cause of the RV variations over nearly 2 
years, in our opinion the covered timespan is not yet sufficiently 
long to corroborate this conclusion and to exclude the outer mas¬ 
sive companion as an alternative explanation. More data are in¬ 
deed necessary over a longer time baseline to cover at least one 
activity cycle of XO-2N and to investigate in greater detail the 
correlation with the RV measurements. 

5.2. The Rossiter-McLaughlin effect 

Dedicated observations of the RML effect for the XO-2N sys¬ 
tem were carried out with HARPS-N during the transit which 
occurred in 2012, December 4-5. The RML observations for 
this system were previously performed by iNarita et~akl (1201 ll) . 
who concluded that the orbit of XO-2Nb is unlikely to be highly 
tilted, but they obtained a poorly determined estimate for the pro¬ 
jected spin-orbit angle, A — 10°±72°. Here we present updated 
results obtained from the analysis of 18 spectra collected at that 



Fig. 10. Phase-folded radial-velocity curve of XO-2N b, after removing 
the quadratic trend. Green diamonds, red squares, and blue circles indi¬ 
cate the HDS, HIRES, and HARPS-N measurements respectively. The 
black solid line shows the best fit model. 
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Fig. 11. Long-term trend found in the RV residuals data of XO-2N after 
removing the orbital signal of XO-2N b, caused either by a long-period 
companion XO-2N c or to the activity cycle of the star. Green diamonds, 
red squares, and blue circles indicate the HDS, HIRES, and HARPS-N 
measurements respectively, while the black solid line displays the best 
fit of the quadratic trend. Residuals of RV measurements with an error 
bar greater than 5 m s -1 , even though considered in the analysis, were 
not included in this plot for a better visual appreciation the long-term 
trend. 


epoch, starting ~ 1 hour before the beginning of the transit up 
to ~1 hour after the last contact. The RML effect was modelled 
with the fitting algor ithm described in lCovino et alJ (120 1 3h and 
lEsnosito et al.1 ( 2014 ). Figure [14] shows the RV measurements 
used for the analysis of the RML effect, folded according to the 
ephemeris calculated in this work, with the best-fit superposed. 
We found a projected spin-orbit angle A = 7°±11°, and derived 
for the projected rotational velocity of the host star Vsin/* = 
1.07+0.09 km/s. The x 1 map in the A - Vsin/* parameter space 
is shown in the same figure. The uncertainties were derived with 
a Monte-Carlo analysis taking into account the uncertainties of 
all the parameters used. Our result provides a stronger evidence 
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Chisquare map 


Fig. 12. Time series of the chromospheric activity index RJj K for XO-2N 
extracted from the HARPS-N spectra. Values of log(RJ |K ) are listed in 
Table [7] We considered only the spectra with a signal-to-noise greater 
than 4 in the 6 ,h order. Superposed is the best-fit linear function, reveal¬ 
ing the presence of a negative trend. 



Fig. 13. Correlation between the RV residuals and the chromospheric 
activity index for XO-2N. Data are extracted from the HARPS-N 
spectra. Superposed is the best-fit linear function, showing the existence 
of a positive trend. 


for a spin-orbit alignment with an uncertainty ~1 times lower. In 
fact, our new estimate of the transit impact parameter partially 
breaks the existing degeneracy between A and the stellar pro¬ 
jected rotational velocity Vsin/*, resulting in a better constraint 
on A. 
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Fig. 14. Upper panel. Radial velocity variations of the XO-2N star in 
correspondence to the transit of the exoplanet XO-2Nb occurred on the 
night between 4 and 5 December, 2012. Data are folded according to 
the transit ephemeris derived in this study [T t .=2,455,565.546480 + 
N*2.6158592 BJD TDB], and superposed is our best-fit RML model 
(red line). Lower panel. The^y 2 map in the A - Vsin/* parameter space. 


orbital inclination angle as measured by transit photometry, and 
A is the projected spin-orbit angle derived from the analysis of 
the RML effect. By calculating sin/* from the projected rota¬ 
tional velocity of the star Vsin/*, the rotation period (41.6+1.1 
days) and the more precise measurement of the stellar radius we 
have derived (0.998+0.033 R 0 ), we obtain sin/* =0.88+0.09. By 
propagating the errors, we finally obtain cos \jj= 0.89+0.12, as ex¬ 
pected very similar to sin/* because the angles i and A are close 
to 90° and 0°. By neglecting the unphysical values for which 
cos this result corresponds to tf/=21^- degrees. Because a 
measurement of the angle i// is difficult to make, out of the 88 
planet host stars for which A was determined, only for nine of 
them if/ was constrained [HI For XO-2N, our estimate is the first 
available. 


5.3. The XO-2N spin-orbit alignment in three-dimensions 

Thanks to the photometric measurement of the XO-2N rota¬ 
tion period, together with the analysis of transit photometry and 
the RML effect, we have the opportunity to provide an esti¬ 
mate of the angle i /r between the stellar rotation axis and the 
normal to the planetary orbital plane. This is given by the for- 
mul a cos i// =cos/*cos/+sin/*sin/cosT (see discussion in Sect.4.4 
of lWinn et al.ll2007l) . where /* is still the inclination between the 
stellar spin axis and the direction of the line of sight, i is the 


6. Stellar activity 

In this work we characterise the stellar activity of the XO-2 com¬ 
panions, both through photometric measurements and spectro¬ 
scopic analysis. The analysis of the photometric time series re¬ 
sulted in an estimate of the spin rotation periods for the two com¬ 
ponents (Sectj4}, we then used the HARPS-N spectra to calcu¬ 
late the chromospheric activity index R[ |K based on the analysis 


http://www.astro.keele.ac.uk/jkt/tepcat/rossiter.html 
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Table 8 . Summary of stellar activity-related information for the XO-2 
system. The values of the activity index log(/?Jj K ) are given as weighted 
mean, and the uncertainty of each parameter is the standard devia¬ 
tion of the mean. The estimates for the ste ll ar rotation per iods corre¬ 
spond to two different relations, those of lNoves et all d!984T ) [N84] and 
iMamaiek & Hillenbrand fc008h [M08], 


Parameter 

XO-2N 
(43 spectra) 

XO-2S 
(63 spectra) 

log (Rind 

-4.91±0.01 

-5.02±0.01 

Exp. stellar rot. period [days] 

39.6±0.6 [N84] 

41.1 ±0.7 [N84] 


41.8±0.8 [M08] 

44.7±0.7 [M08] 

Meas. stellar rot. period [days] 

41.6±1.1 

26.0±0.6 (or 34.5 days) 


of the Can H&K lines. Moreover, an analysis of different CCF 
line shape parameters (asymmetry indicators) was carried out to 
look for periodical modulations and correlations with the radial 
velocities and the R' uw index. 

iFigueira et alJ (12013}) introduced and compared several in¬ 
dicators for the asymmetry of the CCF and shown how they 
can be used to establish whether the radial velocity variation 
of a star is due to line profile distortions induced by, e.g., stel¬ 
lar activity. We have computed the indicators bisector inverse 
span (BIS0, radial velocity difference AV, and velocity asym¬ 
metry V asy for each of the XO-2 components, and a comparison 
between those parameters and the activity index was also per¬ 
formed. We adopted procedure s similar to those of Fi gueira et alJ 
(1 201 3l) and ISantos et al.1 (1201 4l) to compute the asymmetry pa¬ 
rameters with the exception of V asy for which the derivative of 
the CCF with respect to the RV was numerically evaluated by 
means of a Savitzky-Golay filter (e.g., iPress et. all[1992 1. We 
estimated the standard deviation of the indicator AV from the 
standard deviations of the best fit RVs as obtained with the 
Levenberg-Marquardt IDL procedure MPFIT. PRO applied to fit 
a Gaussian and a bi-Gaussian profiles to the CCF, respectively. 
The standard deviation of V asy was computed with a bootstrap 
procedure by taking into account both the photon shot noise and 
the correlated noise of the CCF. 

6.1. XO-2N 

The analysis of the index time series has been presented in 
Sect l5.ll where we discussed the correlation found with the RV 
residuals, possibly indicating that is outlining the activity 
cycle of the star. 

Here we investigate the correlations between the asymmetry 
indicators and the residuals of the radial velocities (i.e. remov¬ 
ing the one-planet fit and the quadratic term), and the results are 
shown in the left panel of Fig. [15J while the right panel shows 
the correlations between the same indicators and the activity in¬ 
dex R' uk - We note that between the CCF asymmetry indicator 
AV and the RV residuals there is a moderate and low-FAP cor¬ 
relation, as quantified by the Spearman’s rank correlation coef¬ 
ficients indicated on the top of each plot. More significant cor¬ 
relations are observed when considering the index instead 
of the RV residuals, in particular for the bisector BIS and AV. 
These seem to indicate that the activity of the star does influence 
the spectral line profile at some extent. 


12 We have verified that between the line-bisector calculated by the 
DRS and the BIS indicator calculated with the method of lFigueira et all 
d2013l ) there is a very high correlation (R=0.99) for both stars, then we 
decided to adopt the latter in our discussion. 
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Fig. 16. Time series of the chromospheric activity index for the star 
XO-2S. Only values of derived from spectra with S/N > 3.9 in the 
6' ,! order were considered. 
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Fig. 18. Can K (upper panel) and H (lower panel) lines of the coadded 
HARPS-N spectra of XO-2N and XO-2S. It is clearly visible a greater 
emission flux in the line cores of the North component, which is indica¬ 
tive of a greater activity of this star. 
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Table 10. Results of the analysis carried out to asses the significance 
of the Spearman’s rank correlations coefficients for the pairs of XO-2S 
datasets shown in Fig. mi 


Method 

R^ k vs BIS 

R^ k vs AV 

R HK vj V o« 

Null hypoth. 

z-score = 0.45 

z-score = 0.37 

z-score = 0.69 


p-value = 32.6% 

p-value = 47.2% 

p-value = 24.5% 

Composite MC 

p = -0.046±0.150 

p = -0.008±0.145 

p = -0.047±0.145 

Original pSpearman 

-0.068 

-0.054 

-0.101 


To better characterize the significance of these correlations, 
we performed the same analysis as described in Section B.l.ll 
and based on two different tests. The results are summarized in 
Table[9] They show that the null hypothesis can be discarded for 
substantially all the pairs, assuming ~1% as the threshold for the 
p-value below which a correlation may be assumed significant. 
The analysis with the composite method suggests that three most 
significant correlations found by simply computing the Spear¬ 
man’s rank coefficient with the actual datasets, i.e. those between 
the RV residuals and the BIS/AV asymmetry indicators of the 
CCF, are likely less pronounced, given the limited number of 
measurements (N=25). Nonetheless, our analysis suggests that 
effects of the stellar activity are influencing the spectral line pro¬ 
files. 

We finally analysed the time series of the CCF asymmetry in¬ 
dicators with the GLS algorithm, to search for signals related to 
the stellar rotation period of ~41 days measured via photometry. 
No significant peaks were found, and in particular no periodicity 
related to the star rotation appears for R^k- 


6.2. XO-2S 


By looking at the average value of the activity index WJ |K (Table 
[8j, this star appears less active than its companion. In Fig. [16] 
the time series of is shown. The internal errors are quite 
comparable to the data scatter, with no clear evidence of a long¬ 
term trend. 

In the same way as for XO-2N, we investigated the relations 
between the asymmetry indicators and the index, and the 
results are shown in Fig. O No significant positive or nega¬ 
tive correlation is found, supporting the scenario of a very quiet 
star. This outcome is confirmed by the results, shown in Table 
ITOl of the analysis performed to asses the significance of the 
Spearman’s rank correlation coefficients, as done for XO-2N. 
For completeness, we also quantified the correlation between the 
index and the RV residuals (after the subtraction of the two- 
planet Keplerian solution) calculated by iDesidera et al . (2 0141 ). 
obtaining a Spearman’s rank coefficient p-0.21. By testing the 
null hypothesis, we found a z-score = -0.94 that corresponds to 
a significant probability p=17.4% that the two datasets are not 
correlated. This result is confirmed by applying the composite 
method, from which we obtain p = 0.20 ± 0.15, indicating lack 
of correlation. 

We also used the GLS algorithm to search for periodicities 
in the CCF asymmetry indexes and /(| |K time series, but no sig¬ 
nificant peaks were detected. For instance, for the BIS and AV 
parameters we found peaks close to half the tentative stellar ro¬ 
tation period derived from photometry (P~26 days). A bootstrap 
analysis resulted in FAP of 0.12 and 0.03, suggesting the the sig¬ 
nals are spurious. No peaks are found close to the orbital periods 
of the two planets orbiting XO-2S. 


6.3. Comparison of the activity levels for the XO-2 stars 


In Table [8] we summarize our activity-related measurements for 
the XO-2 components. We can conclude that, at least in the 
timespan covered by our observations, the star XO-2N appears 
to be slightly more active than the companion. This is also di¬ 
rectly visible when comparing the Ca n H&K lines of the coad¬ 
ded HARPS-N spectra (Fig. HD, which show an emission flux in 
the line cores greater for the North component. Also, the differ¬ 
ences seen in the existing correlations between the CCF asym¬ 
metry indicators and the activity index support the scenario 
that the activity level of the North component is higher. This ap¬ 
pears in agreement with the higher levels of magnetic activity 
in stars with hot Jupiters with respect to their wide binary com- 
panio ns, as found in particular for HD 189733 and Corot-2 (see . 


e.g., Lanzall20Tlt iHusnoo et al.ll2012t IPonnenha eger et alJl2013t 

Popn enhaeger & Wolkll2014l) . Then, the different activity levels 
of the two stars might be due to any influence of the hot Jupiter 
on the activity of XO-2N or to the fact that we simply observed 
the two stars at different epochs of their activity cycles. 


In Table [8] we also show the expected values of the 
spin rotation periods derived from the measurements log(WJ IK ) 
index according to the works of iNoves et~akl (119841) and 
iMamaiek & Hillenbrandl ( j2008i) . It can be seen that the estimate 
of the rotation period of XO-2N derived from photometry is in 
very good agreement with the expected values. 


7. Considerations about the long-term stability and 
evolution of the XO-2 planetary systems 

In a recent work, IWang et ahl (12014l) determined the planet oc¬ 
currence rate in multiple-star systems by measuring the stellar 
multiplicity in a sample of Kepler planet host stars, and com¬ 
paring the stellar multiplicity rate between these and field stars 
in the solar neighbourhood. Interestingly, they find that planet 
formation appears to be significantly suppressed in multiple-star 
systems when the separation between the planet host star and the 
other stellar companion is smaller than 1500 AU, and this dis¬ 
tance could be thus assumed as the radius of influence of stellar 
companions. Nonetheless, the estimate of the planet occurrence 
rate in multiple-star systems is presently a subject of debate. 
iHorch et~aD (12014l) investigated the binary fraction of Kepler ex¬ 
oplanet host stars by observing more than 600 objects of interest, 
and conclude that this appears similar to the rate observed for 
field stars (40-50%), pointing out that the impact of binarity on 
the planet occurrence should be lower than that emerging from 
the work of lWang et al.l d20 14h . 

The analysis of the stellar properties of the XO-2 binary 
points to the scenario that the components have a common origin 
and formed a bound system since their birth. Although the pa¬ 
rameters of the binary orbit are not known, we have performed 
a basic analysis of th e long-term stability of the whole system 
based on the work of iHolman & Wiegerj (11999b . extending to 
the X O-2S system similar considerations don e by iBurke et aD 
(120071) for XO-2N. IHolman & Wiegertl d 19991) modelled the in¬ 
fluence of a stellar companion on the stability of a ’test particle’ 
orbiting a star with an empirical relation (see Eq. (1) therein) 
among the stellar mass ratio /r=m 2 /(mi+m 2 ) (the subscript 1 
identifies the star orbited by a planet), the binary eccentricity 
e, and the semi-major axis of the binary. This relation provides 
an estimate of the critical semi-major axis a c , i.e. the maximum 
orbital distance of a planet to survive after 10,000 binary orbits 
(corresponding to more than 10 9 yr in the case of XO-2). 
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Table 9. Results of the analysis carried out to asses the significance of the Spearman’s rank correlations coefficients for the pairs of XO-2N datasets 
shown in Fig. ED 


Method 

RV resid. vs BIS 

RV resid. vs A V 

RV resid. vs N asp 

RJj K vs BIS 

RJj K vs AV 

R hk vs V^ 

Null hypothesis 

z-score = -2.14 

z-score = -2.93 

z-score = +2.36 

z-score = -4.04 

z-score = -4.32 

z-score = +3.31 


p-value = 1.6% 

p-value = 0.9% 

p-value = 0.17% 

p-value = 0.003% 

p-value = 0.05% 

p-value < 0.002% 

Composite Monte-Carlo 

p - +0.32±0.19 

p = +0.34+0.18 

p =-0.43+0.19 

p = +0.44+0.14 

p =+0.39+0.15 

p = -0.40+0.15 

Original p Spearman 

+0.685 

+0.731 

-0.562 

+0.439 

+0.595 

-0.478 


The work of iHolman & Wieaertl ( 1999 1) considers planets 
that orbit only one component of the stellar system, on initially 
circular and prograde orbits in the plane of the binary. We can 
apply their relation to each XO-2 stars separately, obtaining sim¬ 
ilar results as the mass ratio is p ~0.5 in both cases. By assuming 
the measured projected binary separation a =4600 AU as an es¬ 
timate of the actual semi-major axis of the orbit, and with the 
eccentricity of the binary orbit left undetermined, we can esti¬ 
mate the critical distance in the range of validity of the empirical 
formula, 0<e<0.8. In the case of a single planet orbiting one of 
the components the dependence of the critical semi-major axis a c 
from the eccentricity is shown in Fig. [T9]by the thick curve . Fol¬ 
lowing the discussion in Sect.5 of lHolman & Wiegerli (;1 999 ). we 
can assume half of the a c values found in the case of the three- 
body problem as a rough estimate of the critical distance when 
one star hosts more planets, which is here the case of -at least- 
the star XO-2S. This result is shown by the thin line in Fig.fl9l 
According to it, if the binary has the highest eccentricity consid¬ 
ered here (e=0.8), a planetary companion around XO-2S would 
be in a stable orbit if located at a distance <90 AU from the star, 
while for XO-2N the upper limit is ~ 170 AU. 

A final consideration should be made about the possi ble 
role of the Lidov-Kozai (L-K) mechanism (e.g. Ilnnanen et al.1 
119971: iTakeda et al.ll2008l) in particular on the migration of the 
planet XO-2Nb so close to its parent star. Because the orbital 
parameters of the binary are not known, at this stage we can 
not reconstruct the evolutionary history of this quite old sys¬ 
tem. In particular, we do not have any direct measurement of the 
mutual inclination between them and the binary orbital plane, 
from which the effects of L-K mechanism strongly depends. 
iNeveu-VanMalle et al.1 ( 2014 1 discussed that the L-K mechanism 
likely played a role in shaping the orbits of the two hot Jupiter 
systems discovered around the components of the WASP-94 bi¬ 
nary, this scenario being reinforced by the high projected spin- 
orbit angle 4=151 °±20° and the retrograde orbit of WASP-94Ab. 
Any discussion about the evolution of the XO-2 planetary sys¬ 
tems still remains mostly at the level of speculation. However, in 
our opinion by assuming that the orbital configuration of the bi¬ 
nary did not change greatly with time, the L-K mechanism alone 
can not explain why the XO-2Nb migrated close to its parent star 
and became a hot Jupiter, while XO-2Sb and XO-2Sc are farther 
out, maybe being prevented to migrate because of mutual secular 
perturbations. In fact, the large separation of the two stars leads 
to very long periods of the L-K cycle, well beyond the present 
age of the system. Nonetheless, if the system originated closer 
to the galactic bulge, as discussed in Sect l3.3l this could have 
produced a strong modification of its orbit with time, making 
the study of the perturbing mechanisms very complicated. De¬ 
spite our present inability in assessing the role of the L-K mech¬ 
anism for XO-2 on firm bas is, nonetheless it is inte resting to 
mention the relevant work of lJensen & Akesoiil d2014l) about the 
determination of the protoplanetary disk geometry in a binary 



Binary eccentricity 

Fig. 19. Critical semi-major axis a c for each star in the XO-2 binary sys¬ 
tem, defined as the maximum distance that a single planet, on initially 
circular orbits, could have from its parent star to survive after 10,000 
binary orbits (thick line). The thin line represents a simple estimate of 
a c when one of the s tars hosts more than one p lanet. These results are 
based on the work of lHolman & Wieger ti. ( 1999 ). 

system. Thanks to observations with the Atacama Large Mil¬ 
limeter Array (ALMA), they obtained the clearest picture ever 
seen of protoplanetary discs in a young binary (age < 4 Myr), 
located roughly at the same distance of the XO-2 system. They 
provided evidence for a relative inclination of the discs in the 
range 60°-68°, revealing that both or at least one of them must 
be misaligned with the binary orbital plane, and probably form¬ 
ing an angle far greater than the critical value of ~40° for which 
the L-K mechanism is e xpected to be effective. The results of 
iJensen & Akesonl (120141) suggest that the formation processes of 
a binary can favour those necessary conditions to perturb at some 
extent the planetary orbits since the first stages of planet forma¬ 
tion. 

8. Conclusions 

We have investigated in detail the physical properties of the XO- 
2 stellar and planetary systems. This is the first confirmed case 
of a wide binary whose components, with a T e ff similar to our 
Sun but older and super metal-rich, host planets. This circum¬ 
stance motivated a first comprehensive characterization study. 
Our analysis is based on high-resolution spectroscopic data col¬ 
lected with the HARPS-N spectrograph, and new photometric 
data from several facilities. The main results described in this 
paper can be summarized as follows: 

1. Through a spectroscopic analysis carried out with both the 
equivalent widths and spectral synthesis methods, and the 
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use of stellar evolutionary tracks, we derived new estimates 
of the physical parameters of the XO-2 stars. By performing 
a differential spectral analysis, we measured a positive differ¬ 
ence in iron abundance of XO-2N with respect to XO-2S at 
more than 3<x. When considering this finding in conjunction 
with the existence of a hot Jupiter around the star with the 
higher iron abundance, a very interesting implication is that 
the evolution of this planetary system, assuming the migra¬ 
tion of the giant planet, could have favoured the ingestion by 
the host star of rocky material in the planetary disk, which 
can be revealed as an excess of iron abundance. 

2. Fourteen new transit light curves were collected of the hot 
Jupiter XO-2Nb. From their combined analysis we derived 
updated and very accurate transit parameters. By analysing 
the out-of-transit light curve of XO-2N collected during a 
timespan of ~4 months, we measured for the first time its 
rotation period. We found that the light curve can be folded 
almost perfectly with a sine function of period 41.4+1,1 days 
and that our estimate is well in agreement with predictions. 
We also provide an estimate for the rotation period of XO- 
2S, but our result is inconclusive because the interpretation 
of the periodogram is more complicated than in the case of 
XO-2N. Indeed, the highest peak corresponds to a rotation 
period of 26 days, which is ~ 15 days shorter than that of XO- 
2N and far from being in agreement with gyrochronology 
predictions. However, another peak with almost comparable 
power is found at 34.5 days and, because of aliasing effects, 
for the moment the true periodicity can not be determined. 

3. We analysed the radial velocity time series of the XO-2N 
star, improving the estimates of the orbital parameters of 
XO-2Nb. As a major result, we confirmed the existence of 
an acceleration previously reported in the literature and re¬ 
vealed a turn-over in the trend. We discussed the implications 
of this finding in terms of the presence of an outer compan¬ 
ion for which we determined a lower limit for its mass and 
orbital period. Nonetheless, we found a significant positive 
correlation between the residuals of the RV and the activity 
index which would imply that the observed acceleration 
is rather caused by the intrinsic activity of the star and related 
to its solar-like activity cycle. However, despite we found 
relevant and convincing indications for a stellar activity cy¬ 
cle as the real cause of the observed long-term RV variation, 
more observations would be advisable to definitively clarify 
the origin of the RV trend. Finally, by combining the results 
from the analysis of the transit light curves and of the radial 
velocity time series (archival data and new HARPS-N mea¬ 
surements), we refined the main physical parameters of the 
planet XO-2Nb, in particular the mass, radius and density. 

4. We modelled the Rossiter-McLaughlin effect for XO-2N, 
providing a better constrain on the projected spin-orbit angle 
A-l°±\\°, according to which the system is likely aligned. 
By combining the result on A with that on the stellar rotation 
period, we could derive for the first time an upper limit of 
39° for the true angle t ft formed by the stellar rotation axis 
and the normal to the planetary orbital plane. 

5. We investigated the stellar activity of the XO-2 components 
by analysing a set of asymmetry indicators of the spectral 
cross-correlation function and the chromospheric activity in¬ 
dex We conclude that the North component was more 
active during the timespan of our observations, and the dif¬ 
ferent activity levels might be due to any influence of the hot 
Jupiter XO-2Nb on its host star or to different phases of their 
activity cycles. 


The XO-2 system should be still regarded with interest, in 
particular for the evidence of long-term trends in the RVs of both 
the stars for which a clear explanation is still not possible. Should 
these signals be caused by additional companions in outer orbits, 
this would make the system still more interesting to be studied as 
a peculiar laboratory for testing the theories of planet formation 
and evolution. We are collecting new data to improve the char¬ 
acterization of the whole system and provide insights to several 
open questions. 
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X02-N 



X02-S 



Fig. 2. Iron abundance of the XO-2 stellar binary derived through the analysis based on equivalent widths as a function of the excitation potential 
and reduced EW. Empty and filled circles represent the Fe i and Fe n lines, respectively. Solid lines show the location of mean log n {Fe i), while 
dashed lines represent the standard deviation ±lcr. The dotted line is intended to guide the eye to the solar iron abundance at log/;(Fei) 0 = 7.53. 
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Fig. 4. Magnitude differences between the star XO-2N and its companion XO-2S as measured in BVRI passbands from data of the TASTE project, 
APACHE survey, and Serra La Nave Observatory. These differences were evaluated from out-of-transit photometry. 
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Fig. 8. Transit light curves of the planet XO-2Nb collected with the Copemico (182-cm) and Schmidt 67/91-cm telescopes (Asiago Observatory) 
within the framework of the TASTE Project (first ten plots and the last one) and, independently, at the IAC-80 telescope (Teide Observatory, Canary 
Islands; remaining plots). All these observations were carried out in the Cousin/Bessell R-band. The light curves are plotted both in their original 
cadence (left diagrams), and binned over 2-rnin intervals (right diagrams). For each transit, the RMS of the residuals from the best-fit model is also 
tabulated. A composite light curve has been obtained from all the individual time series and used for determining the transit parameters shown in 
Table[5] as discussed in the text. 
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Fig. 15. Correlations between the asymmetry indicators and the residuals of the radial velocities (left panel) for the star XO-2N, and between the 
same indicators and the activity index R' UK (right panel). We considered only values of R' ]]K extracted from spectra with a signal-to-noise greater 
than 4 in the 6 th order. 
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Fig. 17. Correlations between the asymmetry indicators and the chromospheric activity index for the star XO-2S. Only values of R' HK derived 
from spectra with S/N > 3.9 in the 6 ,; ' order were considered for the analysis. 


Article number, page 26 of 1261 

























































































































































